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This thesis is an application Df parameter sensitivity 
analysis to aerodynamic parameters of a Bank-to-Turn 
missile. In the development a brief review of trajectory 
sensitivity theory is presented. A linear analysis is 
performed using an Uncoupled Pitch Channel Autopilot and a 
Coupled Roll-Yaw channel Autopilot of the missile taken as 
model. Finally, a nonlinear analysis is given for the system. 
Comparisons between the linear and nonlinear cases are 
outlined. 
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I- p arameters sensit ivity a nalysi s op a BANK-TO-rORN MISSILE 

The determination of changes in system performance due 
to changes in parameters is of great importance in engi- 
neering analysis and design. Sensitivity guest ions arise 
when model uncertainty is present, or when a range of oper- 
ating conditions is contemplated. 

The guestions of parameters sensitivity particularly arise 
in the fields of engineering where models are idealized, 
inexactly identified, or the systems themselves are subject 
to unpredictable changes with tiae due to en viromental, 
material property or operational influences so that there is 
always a discrepancy between the physical reality and the 
mathematical model. 

Sensitivity analysis provides the engineer with methods for 
investigating or minimizing such parameter deviations. 

In general, the diagram of a system can be represented 
by a single block as given on Fig. 1.1. 











SYSTEM 









Figure 1.1 General Representation of a Dynamic System. 

From a mathematical point of view, what we call a 
system is the explicity or implicity given relationship 
between the input signal u and output signal y. In general, 
u and y can be vectors. The character of this relationship 
is commonly called the structure of the system. 
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For example, the structure of the system may be 
characterized by the order of a differential or difference 
equation, linearity or nonlinearity, the order of the numer- 
ator and denominator of a rational transfer function and the 
quantitative properties of the system parameters. 

Typical parameters are initial conditions, time- 
invariant or time-variant coefficients, natural frequencies 
and samplinq periods. 

The chanqe of the state or the change of the output 
variable with time, can be caused by; (1) the influence of 
input signals, (2) the change of parameters. These quanti- 
ties are shown in Fig. 1.2. 




Figure 1.2 Quantities Affecting the Dynamics of a System* 

This work addresses the application of the sensitivity 
analysis to the variation of aerodynamic parameters of a 
Bank-to-Turn missile. 

PARAHETER SENSITIVITY is the affect of parameter changes on 
the dynamics of a system, say, the time response, the state, 
or any other quantity characterizing the system dynamics. 

In order to have a realistic model to work with a NASA 
Contractor Report [Ref. 2] was adopted as a reference for 
application of the sensitivity analysis on the variation of 
aerodynamic parameters of a Bank-to-Turn missile. Appendix 
A, B and C give the detailed aerodynamic assumption of the 
system. 
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One brief explanation of the sensitivity theory and 
analysis is given as a guideline to better uniersrand the 
subsequent work. 

An analysis of the sensitivity of the aerodynamic parame- 
ters applied to the case of the linear uncoupled pitch and 
coupled roll-yaw autopilot is performed. 

The pitch and roll- yaw autopilots are coupled as given 
in Fig.C.1 in the appendix C. A nonlinear parameter sensi- 
tivity analysis is performed. 

Conclusions about the influence of parameters of concern 
are given and a comparison between the linear and nonlinear 
case is done. Comments and recommendations for the future 
are outlined in order to delineate the continuity of the 
present work. 



22 



I 

i 

I 



J 

n 

I 



II. SENS ITIVIT Y THEORY 



A. INTRODOCTIOH 

The basis of all sensitivity considerations in the case 
of time-invariant parameter variations is the so-called 
sensitivity functions. Dynamic systems can be characterized 
in several ways: in the time domain, in the frequency domain, 
or in terms of a perfomance index. There is evidently an 
adequate number of ways to define the sensitivity functions 
of a dynamic system. The definition that is actually used 
depends on the form of the mathematical model as well as on 
the purpose under consideration. 

The sensitivity functions can be classified into the 
following three categories: 

(1) Sensitivity functions in the time domain 

(2) Sensitivity functions in the frequency domain, and 

(3) Per fomance-index sensitivity 

In this chapter we will outline just the sensitivity func- 
tions in the time-domain of the continuous systems. 
Aplications of this analysis follow throughout this work. 

The reader can find detailed informations about the 
others two categories of sensitivity functions in the 
[Ref. 1], 

B. TRAJECTORY SENSITIVITY FUNCTION OF CONTINUOUS SYSTEMS 

A continuous, possibly nonlinear system of nth order 
can, in general, be described in the state space by a vector 
differential equation as seen in Eqa.2.1. 
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X = f(X,t,U,'<o) 



( 2 . 1 ) 



, X(t^) = XO 

Here X is an n x 1 state vector, f an n x 1 vector func- 
tion, "u an input vector, a nominal r x 1 parameter 
vector, and X is the n x 1 initial condition vector or 
initial state. Eqn.2.1 is called the NOMINAL STATE 
EQUATION. 

Assuming that the parameter vector deviates from the 
nominal value by we have the Eqn.2. 2. 



X 



f(X,t,U, OC), X(to) 



XO 



( 2 . 2 ) 



with the initial conditions X® unchanged. This equation 
is called the ACTUAL STATE EQUATION. 

Now it is assumed that Eqn.2. 2 has a unique solution, X= 
X(t,^) for all admissible initial conditions and parameter 
values. 

X is of course a function of u , X° and t as well. 
However , this dependence is not needed for the following 
considerations and will, therefore, be dropped for ease of 
notation. Furthermore, the solution X is assumed to be a 
bounded continuous function in t and ^ . 

If the parameter takes on its NOMINAL value <x-o » the 
nominal solution X^ = X(t,o^) is obtained. If, on the 
other hand the ACTUAL solution is given by X ”= X(t, ) , 
then the parameter- induced change of the state vector is 
given by 



X(t,<K) = X(t, X) - X(t, CK^) 



(2.3) 
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k first-order approximation of ^ X can be obtained by 
using Taylor expansion in the form of Eqn.2.4. 



A X (t , o< ) 



2 

Jr1 



(2.4) 



This equation can be viewed as a definition of the 
parameter-induced trajectory deviation. 

The subscript oc© shall indicate that the partial deriv- 
ative expressed by ”D is taken at nominal parameter values. 

Let the state X of a continuous system be a continuous 
function of a time- invariant parameter vector o< = 

{ o<, , c <2 • Then the partial derivative can be 

defined as 

- - 

A;(t,^«) = VS”*”” 

f 2^ . . . / n 

j— 1 f 2g • • • f T 

fiqn.2-5 is called the trajectory sensitivity vector with 
respect to the jth parameter. 

Note that the trajectory sensitivity vector is of the 
same dimension as the state vector, namely, n. Its compo- 
nents are the TRAJECTORY SENSITIVITY FaNCTIONS as 



T .. . A (t,ot) 

"J ^ 




( 2 . 6 ) 



Eqn.2.6 is the partial derivative of the ith state vari- 
able in relation to the jth parameter as 
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(2.7) 






.X 









— _j 

D «-o 



Hence, all n x r trajectory sensitivity functions form 
the trajectory sensitivity matrix as given in Eqn.2.8 or 

2.9. 



X 





A 'll X 




r "t) "1) xi 



A 



0 



^Xn ^ ^ Xn 



( 2 . 8 ) 



(2.9) 



The columns of A are the trajectory sensitivity 
vectors . Here X is the Jacobian matrix of the state 
vector with respect to the parameter vector «. , taken at 

the nominal parameter values. 

With these definitions the parameter-induced change of 
the trajectory can be taken as 



ir __ 

Zilct,^) = A (t, <x^) = X (2.10) 



Where o< = , which is the iCTOiL parameter 

vector of the system. 
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C. TBAJECTORY SENSITIVITY EQUATIONS OF CONTINOOOS SYSTEMS 

Lets consider the general continuous system described as 
previously by the state equation (Eqn.2.1). 



X = f (X,t,U, X), X(t^) = xo 



( 2 . 11 ) 



Where X denotes the n- dimensional state vector, ^ the 
r-dimensional parameter vector, f an n-dimensional vector 
function, and u the input vector independent of c><. . it is 

assumed that the continuity conditions are fulfilled and 
that o<. is time-invariant. 

Considering c< -parameters and taking the partial 
derivative of X (Eqn.2.11) with respect to «<j' , one obtains, 
by the application of the chain rule 



~i) X 'll X X® 



The derivative of the initial conditions vector X° with 
respect to c<j* is zero, since X° does not depend on X . 

If X is r-dimensional, there are r equations of the 

form of Eqn.2.1 2. If we now interchange the sequence of 
taking the derivative wirh respect to time t and®*<’j’, and 
then let X approach to ^ one obtains the following 

equation 




1.1 
^ X 




X'o)= 0 

^ «» 



(2. 13) 
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Here is the trajectory sensitivity vector 

-J Q 

with respect to the jth parameter. 

The Eqn.2. 13 is called the sensitivity equation in the 
state space or the TRAJECTORY SENSITIVITY EQOATION. 

The above shows that for c< -parameters all initial 
conditions of the trajectory sensitivity equations are equal 
to zero. 

Now consider the output vector equation, as given by 



Y = g(X,t,u,<X.) 



(2. 14) 



In a procedure similar to that above one obtains the 
algebraic sensitivity equation as seen in Eqn.2. 15. 


















i;. . 



-• 












(2.15) 



the output sensitivity 



vector = 
A.- 



_ _ which relates 

0 y 

to the trajectory sensitivity vector ^j' . This 
equation'^is called the VECTOR OUTPUT SENSITIVITY EQUATION. 

Using the trajectory sensitivity matrix A and the 
output sensitivity matrix , the above result can be 

rewritten in the following general form 



^ 'T)4 T 

A = — r. — U A + 









1)4 T 

1 _ .X» = 

OoC ®<o 



(2. 16) 



cr = X ♦ 

1>X- 



^ OL °<0 



(2. 17) 
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These equations are called the STATE SENSITIVITY 
EQUATIONS of the system. It is seen that these equations 
are linear whether the original system is linear or nonlin- 
ear. 

If, in particular, the original system is Linear, the 
state equations, take the form 



X = AX + Bu ,X(to) =XO 



(2. 18) 



y = c X + D u 



(2.19) 



Where, in general, 

r=T(o<.) ,'b = ¥(o<.) ,c = c(o<) ,d = ¥('^), 

T = F(t,c< ) , amd 7 = 7(t,6< ) . 

Note, however, that u is not a function of «<. if u is 
defined as an external input of the system. 

Now taking the partial derivatives with respect to , 

reversing the order of differentiations with respect to 
time and o4j , and letting ^ approach to , the 

TRAJECTORY SENSITIVITY EQUATIONS are obtained 



X; . s, X\ * ^ 1 x« * 









I- u(t) , ACl:o) = 0 (2.20) 






Where A^ = k{c<o) , X^, = X(t, o^q) , and j= 1, 2,..., r. 

The initial condition vector X -Lt:©), is again zero since 

sj 

X(to) does not depend on c< . 
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By similar procedure the VECTDR SENSITIVITI EQUATION 
becomes 

~ -- 'D C - 'I) ^ - 

LXo ♦ l_u(t) (2.21) 

® lio<^ oc^ 

Where = C ( ck^ ) and j=1,2,..., r. 

In the case of linear systems, the vector sensitivity 
equations have the same A matrix as the nominal state equa- 
tions and hence the same characteristic matrix si - A. 

They differ from the nominal original state equation only 
in the driving function and the initial conditions. The 
latter are all zero. The driving functions can be obtained 
by solving the nominal state equations. 

Here, in the linear case, the sensitivity equations can 
be joined to the system equations , forming the so-called 
COHBINED SYSTEM as given in Eqns.2.22 and 2.23. 







> 

a 

0 

i 








1 

CPJ 

0 




• 




"I> A 1 ^ 






+ 


■^8 




. X . 














*<0-* 



u.c-t) 



(2. 22) 







■ Co O 




Xo 




Do 




* 


— 


1 






+ 


'Do 








‘"'’J 
















(2.23) 



The simultaneous solution of these differential equa- 
tions using the same standard procedure for each of the 
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above matrix equations yields the output vector, the state 
vector, and the corresponding output and trajectory sensi- 
tivity vectors , and Aj • If there are r parameter 

variations, r systems of equations of the above type have to 
be solved. Since the homogeneous part of the original 
differential equation is identical with the homogeneous part 
of the sensitivity equations with respect to all parameters. 

A graphical interpretation of Eqns.2.22 and 2.23 is 
given in Fig.2.1. 

Since the driving function of the sensitivity model 
contains the nominal state, the measuring circuit for the 
trajectory sensitivity functions consists of a connection of 
the nominal original system and the sensitivity model as 
illustrated by Fig.2.1. 

If the actual input u is applied to this structure the 
trajectory and output sensitivity vectors Aj and <l^- can be 
measured at the points 1 and 2 respectively in the Fig.2.1. 

In order to measure all r sensitivity vectors simultane- 
ously, r sensitivity models are required. 

D. STEOCTURAL METHOD 

The basis for the measurement of the trajectory sensi- 
tivity functions is the STRUCTOBAL INTERPRETATION Of the 
trajectory sensitivity equation. 

The physical system represented by the trajectory 
sensitivity equations is called the trajectory sensitivity 
model or the sensitivity model in the stare space. 

Regardless of the nature of the original system the 
sensitivity model in the state space is always linear. The 
graphical illustration of sensitivity equation can be given 
by Fig. 2. 2. The system described by the sensitivity equation 
is refered to as the sensitivity model of the original 
system. For each output, a system has as many sensitivity 
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models as parameters of interest. Both the nominal original 
system and the corresponding sensitivity models form the 
COMBINED SYSTEM. The sensitivity model is always linear. If 
Y(t,^) is the output of a syste^ the corresponding sensi- 
tivity functions (t,c<p) = outputs of 

the corresponding sensitivity models. Consequently, in order 
to measure the output sensitivity functions simultane- 
ously, the nominal original system and the sensitivity models 
have to be measured at the outputs of the sensitivity 
models. 

The double frame used for the original system in Fig. 2.2 
is to indicate that, in general, the system may be nonlin- 
ear, whereas the sensitivity models are always linear. In 
the nonlinear systems the Eqn. 2.12 applies and different 
programs have to be set up for the original and the sensi- 
tivity equations. 

For some applications, the sensitivity functions with 
respect to all parameters are required simultaneously. If 
there are r parameters of interest, r sensitivity models 
would be needed. This implies a rather extensive computer 
time . 

A method of determining all output sensitivity functions 
of a system simultaneously by a single sensitivity model is 
available , which is also called the method of SENSITIVITY 
POINTS. This method has application just for linear systems. 
Detailed explanation about the application of the sensi- 
tivity points theory can be found in the [Ref. 1]. 

E. OOTPOT SENSITIVITY PONCTION OP CONTINOOOS SYSTEMS 

Consider the input-output behavior of a continuous, 
possibly nonlinear, single-variable system described by an 
ordinary differential equation of the type 
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} — 0 



( 2 . 24 ) 



f{ y 



in) (n-<; 



with the initial condition y(tj,) = y. o, i=0, 1,...n-1. y 
denotes the output signal, t the time, and cx a single 
time-invariant or slowly varying -parameter that has 
nominal value • 

In general, f is a function of the input u as well. 
However, if u is an external input which does not depend 
on o<^ , the dependence of f on u is not relevant in further 

considerations , 

Let us suppose that the above NDMINAL differential equa- 
tion has the unique solution given below 



lo = Y(t, ®<o) 



( 2 . 25 ) 



Which one shall call the NOMINAL SOLOTION. 

Let us now assume that the parameter changes 

from to <K= where ^ is time-invariant or 

slowly varying with time. <X. is called the actual parameter 
value. The corresponding ACTUAL DIFFERENTIAL EQUATION can, 
then, be written as 

in) 

f — {iJ,y ,...,y,t,o(.)=:0 (2. 26 ) 



Note that by this change of o<o into oc. the initial 
conditions remain unchanged, namely y(t^) = y.®* 
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The corresponding solution is given as 



y = y (t, o< ) 



(2. 27) 



Which we shall call the ACTOAL (or PERTURBED) SOLUTION. 

It is assumed that y(t, o<) is of the same type as 
y(t» ) » and y(t, o<) deviates infinitesimally from y(t,°<o) 
if deviates infinitesimally from • The conditions 

for fulfilling this requirement are given in the mathemat- 
ical literature. For our purpose, it is sufficient to know 
that y is continuous in if f is continuous in y which is 

true for all continuous systems and ^ 0. 

With the above assumptions, the actual solution y (t, 

+ ^ae.) can be expanded into a Taylor series around °<o , 
yielding the Eqn. 2.28. 



If << , the Taylor series can be truncated at the 

linear term. 

This gives the Eqn. 2.29. 



For finite values of ^ , this expression can be 

considered a first-order approximation of y(t,o<.). 

The actual output can be written as 



y(t,<x) = y(t,og + 




--fl +. 



(2. 28) 



y(t,of) = y(t,(Xo) + 1 ^ 



(2.29) 
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y(t,<x) = y(t,c<o) ♦ ^ (t,o<o) a 



( 2 , 30 ) 



Where y(t,«^o) is the nominal output and (T ( 
the parameter-induced output error. 

With (T (t, o<.o) defined | . Ttie 

induced output error is in these terms defined as 



y (t,^ ) 



A 



^ (t,«o) 



A oc 



parameter 
( 2 . 31 ) 
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Figure 2.1 Sraphical Interpretatiaa of Eqns.2.22 and 2.23. 
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Figure 2.2 Strutural Diagram of the Combined System. 
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III. APPLI CAriON OF SENSIT IVITY TO UNBAR SY STE MS 

A. IHTRODOCTION 

The sensitivity theory described previously will be 
applied to the case of an uncoupled pitch autopilot and to 
a roll-'yaw coupled autopilot presented in the appendix C. 
For each case, the linear equations of the nominal system 
are presented in state variable form and the correspondent 
trajectory sensitivity equations are, derived. A sensitivity 
analysis of the systems is performed. 

B. ONCOOPLED PITCH AOTOPILOT ANALYSIS 

Linear Equa tions of the Nominal System 

From the block diaqram of the Fig.B.1 in the 
appendix B one can obtain the following nominal state equa- 
tions of the uncoupled pitch autopilot. 



XI 



= C2 A3 X2 + C2 A4 X3 



(3.1) 



X2 = XI - K Cl A1 X2 - K Cl A2 X3 



X3 = - C3 X3 + C3 Conv X6 



X4 = - Cl C4 A1 X2 - Cl C4 A2 X3 -C4 X4 



(3.2) 



(3.3) 



(3.4) 
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X5 = C7 X4 - C5 X5 - C6 NZC 



(3.5) 



X6 = C9/Conv + C2 C8/conv A3 X2 



(3.6) 



♦ C2 C8/Conv A4 X3 - C7 C8 X4 ♦ (C5 C8 - C9) X5 
+ C6 C8 HZC 



For the purpose of using the linear method these 
state equations must be presented in matrix form as given by 







o 






o 


o 


•o 






• 


o- 






1 




-KC^A^ 


o 


o 


o 








o 






o 


o 




o 


o 


CjCOhv 






+ 


o 






o 


“ ^ 




-^4 


o 


o 








o 






o 


o 


o 






o 










# 


i 


C9 


CrCe As 


CfcCe 










^6 






6 


1 


Cony 


Qonv 


conv 













Nzc 



(3.7) 



The correspondence of the state vectors is: 

XI = c| , X2 = , X3 = , X4 = X / X5 = Y , X6 = S 

Definition of the constants Cl through C9 are given 
in appendix B. The parameters of interest for this system 
are given by A1, A2, A3, and A4 which are : 



A1 = C., , A2 = , A3 = C 

This nomenclature is 
sensitivity theory and to 
carrying symbols and constants 



, A4 = . 

used here to easily apply the 
avoid the inconvenience of 
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2. Sensit ivity Equati ons 

To apply the procedure developed in chapter 2, 
consider the trajectory sensitivity equation (Eqa.2.13). 






= A 



X. * |_x 

^ (X • <V . 



+ — |_u(t), 



(3.8) 



From Bqn.3.7 one can see that 



V 



o 






O 


O 


■O 


1 


-KiC,A, 


-KC.A^ 


O 


O 


o 


o 


O 




O 


o 


CjConv 


o 






-^4 


o 


O 


o 


o 


o 


Cf 




O 


ce 

conv 


conv 


C.C9 A. 
conv 






% O 



(3.9) 



E>- [ O O 



- C 



c 



Ce Cs 3 



(3. 10) 



Hi D ft 

The partial derivatives -^-^-.-|^and -l ^^ are 

evaluated considering the parameters of interest A1, A2, A3 
and AU, that are respectively the aerodynamic coefficients 
present in the pitch channel. 
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Applying the partial derivatives with respect to the 
parameters we have the following matrices 



HA 



o 


o 


Ca 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


O ' 


o 


6 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


Ca C© 
Ce»f\v 


o 


o 


o 



( 3 . 11 ) 





o 




o 


o 


o 


o 




o 


o 


o 


o 


o 


o 




o 


o 


o 


o 


o 


o 


■JA . 


o 


o 


o 


o 


o 


o 




o 


o 


o 


o 


o 


o 




o 


CzCq 

conv 


o 


o 


o 


o 



( 3 . 12 ) 



1) Ajj 



o 


o 


o 


o 


o 


o 


o 


o 


-KC, 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 




o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 



( 3 . 13 ) 



"b A 
7) A/ 



o 


o 


o 


o 


o 


o 


o 




o 


o 


o 


o 


9 


o 


o 


o 


o 


o 


o 




o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 



( 3 . 14 ) 
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because B is independent of the parameters of interest. 
In terms of components one can see 







(3. 15) 



i=1,2,...,6 j=1,2,3 and 4 

Here, four models are necessary to study the effect 
of parameter variations. These models are shovrm in Fig.3.1. 

For instance, when one apply the theory for one 
parameter of interest, say A1, the sensitivity equations can 
be obtained from 



" • 








* 
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^1 


?11 








11 




























^2 


^3 
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o> 
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+ 




^3 


a" 
















1 





(3. 16) 



Similar procedure can 
ters A2, A3, A4. 



be done 



for the other parame 




Figure 3.1 Nominal and Sensitivity Models. / 
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C. SENSITIVITY ANALYSIS 



A computer program using [Ref. 4] was writtsn in order 
to simulate simultaneously the nominal and the sensitivity 
models as shown in the appendix D. A step was applied as 
input of the nominal system. 

The number of equations solved are 6 for the nominal 
model and 24 for the sensitivity model. Each parameter was 
varied simultaneously 1056 from all the nominal values. 



D. ANALYSIS OP -PARAMETER VARIATIONS 



The results are plotted in Fig. 3.2 to 3.6 and Table I 
Each state variable output is plotted separately with the 
correspondent four output sensitivity functions. By means of 
these plots, the following observations can be made: 

The plot of /I,; = • j=1»2,3, and 4 in fig. 3.2 

indicates that a parameter change Ajg> ( j= 1 ,2, 3 , and 4) 
primarily affect the rise time and overshoot of X1(t) since 
>i11, /112, /113, A14, are largest at the time where these 



effects in XI occur. ^ 

The plot of / 3 . 2 ^= j=1»2,3, and 4 in Fig. 3. 3 

indicates that X2(t) is little affected due to ^ Aj^ and 
A A 20 and strongly affected in the rise time due to param- 
eter variations ^Ag^j and ^ A^c>< ^^e overshoot and steady 

state are little affected due to parameter variations 



/h A j (j=1,2,3, and 4) . 

^ X 3 

The plot of ' j=‘l#2,3, and 4 in Fig. 3.4 

indicates that X3(t) is little affected in the rise time due 
to parameter variations A (j=1,2,3,and 4). The overshoot 
and steady state are little affected due to parameter varia- 
tions ^A|o and Al .20 and strongly affected due to A A 30 and 



A Aj-O. 
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The plot of j = j=1»2,3, and 4 in Fig. 3.5 

indicates that X4(t) is little affected in the rise time and 
overshoot due to parameter variations ^Aio 
strongly affected due to ^ and 

The plot of Atr‘,= r j=1,2,3, and 4 in Fig.3.6 

5J A_i Aj o 

indicates that X5(t) is little affected in the rise time due 
to parameter variations A A^q (j=1»2,3, and 4). The over- 
shoot is little affected due to A A^^ and ^ A^q and strongly 
affected due to Aw a 3 ^ and a. A^q. 

The plot of X(,j= j=1,2,3, and 4 in Fig. 3. 7 

indicates that X 6 (t) is little affected in the rise time due 
to parameter variations A A j ^ (j=1,2,3, and 4). The over- 



shoot and steady state are little affected due to A A 



(O 



and 



A A 2 o strongly affected due to A a^q and ^ 

Table I shows the above analysis in a condensed way to 
give a general picture of all states and sensitivity func- 
tions with the correspondent effect as a function of time. 



E. PAHAMETEH-IHDOCED OOTPOT ANALYSIS 

If A«< « Xo ,the actual output can be written as 
(Eqn. 2.29) . 



y(t,x) = y(t,o<^) 






(3.17) 



The computer program as shown in the Appendix D gives 
the actual output when 10%, 3015 and 40% of variation from the 
nominal value of each parameter is assumed. Figs. 3.8 
through 3.13 show the actual output when each parameter is 
varied 10 % from the nominal value. 

Fig. 3 . 8 indicates that the overshoot and rise time of 
the actual and nominal output of XI are strongly affected 
which is in agreement with the previous analysis. 
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?ig.3.9 iadicates that the steady state and overshoot of 
the actual and nominal output of X2 are primarily affected 
which is in agreement with the previous analysis. 

Fig. 3. 10 indicates that the staady state and overshoot 
of the actual and nominal output of X3 are primarily 
affected which is in agreement with the previous analysis. 

Fig. 3. 11 indicates that the overshoot of tha actual and 
nominal output of X4 is primarily affected which is in 
agreement with the previous analysis. 

Fig. 3. 12 indicates that overshoot of the actual and 
nominal output of X5 are primarily affected which is in 
agreement with the previous analysis. 

Fig. 3. 13 indicates that the steady state and overshoot 
of the actual and nominal output of X6 are primarily 
affected /having little effect in the rise time which is in 
agreement with the previous analysis. 

Figs. 3. 14 through 3.19 and Figs. 3. 20 through 3.25 show 
respectively the actual output for 30% and 40% of the 

nominal value. 

From the plots one can see that for small parameter 
variations the parameter-induced output error is negligible 
and when the variation becomes larga as 30% or 40% one note 
that the error becomes pronounced and that modeling is 
starting to break down. This agrees with the assumption 
made in the derivation of the Eqn.2.28. 

P. COOPLED ROLL-Ika AUTOPILOT ANALISIS 

1. Linear Equ at ions of the N omi n al Syste m 

From the block diagram of the Fig.B.1 in the 

appendix B one can get the following state equations of the 

coupled roll-yaw autopilot. 
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XI = B conv (A4 X3 + A3 X8 ♦ A5 XII) 



(3. 18) 



X2 = C Conv (A7 X3 + A8 X8 ♦ A6 XII) 



(3 



X3 



XI - ( ALPHA B/Conv) X2 



+ KB A A2 X3 
X4 = - 8 X4 - 



+ KB A A1 XII 
17.6 X12 + 17.6 PHC 



(3 



(3 



X5 = - (0.755/Coqv) X2 - C D A7 X3 (3 

♦ (0.755 - 8 D) X4 - 5 X5 - C D A8 X8 - C D A5 XI 1 
- 17.6 X12 +17.6 PHC 

X6 = C E A7 X3 - 6 X6 ♦ B C A8 X8 + C E A6 XII (3 



9 

X7 = - F KC (0.755/Conv) X2 (3 

- (D + E) F KC C A7 X3 + P KC (0.755 - 8 D) X4 
+ KC(15 - 5 F) X5 + KC(6F -15) X6 - 15 X7 

- (D + E) F KC C A8 X8 - (D + E) F KC C A6 XII 

- F KC D 17.6 X12 > F KC D 17.6 PHC 

X8 = 188.4 Conv X7 - 188.4 X8 (3 



.19) 



. 20 ) 



. 21 ) 



. 22 ) 



. 23) 



.24) 



.25) 
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X9 = (K1 k A2/t1) XI - X9 / T1 + (K1 k A1/T1) XII (3.26) 



X10 =- (K2/Conv) XI - (K2 H ALPHAB/Conv X2 (3.27) 

(k2/10) (K1 A A2/T1) ♦ B A4 - H ALPHAB C A7 X3 
♦ (K2/10) (B A3 - H ALPHAB C A8 X8 + K2(1 

- 1/(10 T1)) X9 + (K2/10) ( (K1 A A1) /T1 + B A5 

- H ALPHAB C A6) XII 

XII = 188.4 Conv X10 - 188.4 XII (3.28) 



X12 = X2 / Conv (3.29) 



Eqn.3. 30 gives the matrix representation of the 
state variables above mentioned. 




(3.30) 
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The corr espondence of the states is; 

X1 = r,X2 = p,X3 = fi , 

X4 = X , X5 = T1 , X6 = XI , 

X7 = ^R.c , X8 = 6«. , X9 = Y , 

X10 = Syc ' X11 = Sy , X12 = 4> 

The definition of the constants Cl throagh C25 are 
given in appendix B, The parameters of interest for this 
system are given by A1 through A8 which are : 



A1 = Cy, , A2 = Cy^ , A3 = , A4 = C„ , A5 

A6 = C- , A7 = C,„ , A8 = Co . 

<Sy 

2. Sensitivit y E quati ons 



c 

n 6 y 
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As showed in the previous analysis here one can 
apply the same procedure using the TRAJECTORY SENSITIVITY 
EQUATION as given in Equ.2. 13. 




o j 



't) A - 
--C----LX0 

Vo 






Uu(t) 

'<0 



(3.31) 



For the purpose of this procedure again one can see 
in Eqn.3.30 the correspondent matri£es A^, and B^. ^ 

The partial derivatives and the are 

evaluated considering the parameter A of interest. In this 
case they are A1 ,A2 , . .. , A8, that are respectively the aero- 
dynamic coefficients present in the coupled roll-yaw auto- 
pilot. 

Once again, applying the partial derivatives with 
respect to the parameters one can obtain eight matrices 
respectively as found , similarly, in the previous case of 
the uncoupled pitch autopilot. 

Here, one can see that xr — r = 0 because B is inde- 

o 

pendent of the parameters of interest. 

In terms of components one obtains 
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(3. 32) 



i=1 ,2, 12 

j=1,2, 8 

Here, one can see that eight models are necessary to 
stady the effect of parameter variations. These models are 
shown in Fig. 3. 27. For one parameter of interest, say A1, 
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(3.33) 



Similar procedure can be male for the others parame- 
ters A2 through A8. 



G. SENSITIVITY ANALYSIS 

In order to simulate simultaneouly the nominal and the 
sensitivity models, a computer program was written as shown 
in the appendix E. 

For analysis purpose, each parameter was varied simulta- 
neously 10% from the nominal value. 
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Analy sis of o< -Parameter Variations 

^ X, 



The plots of A,,’ = — Ia' » j=1,2,.,.,8 in 

-* o /\i o 

Fig. 3.27 ana 3.28 indicate that a parameter change 
(j= 1 , 2 ,. .. , 8) primarily affect the overshoot of K 1 (t) . 

The plots of A -2 |‘ = Ia." » j=1»2,...,8 



in 



1>Aj 

Fig. 3.29 and 3.30 indicate that a parameter change AAJf,, 
(j=1 , 2,. .. , 8) primarily affect the overshoot of X2(t). 

The plots of j=1,2,...,8 in 

Fig. 3.31 and 3.32 indicate that a parameter change 
(j= 1 , 2,. . . , 8) primarily affect the overshoot of X3(t). 

The plots of A^j = j=1,2,...,8 in 

Fig. 3.33 and 3.34 indicate that a parameter change 
(j=1 ,2 ,. . . , 8) primarily affect the overshoot of X4(t). 



The plots of Agj = j=1r2,...,8 in 

Fig. 3.35 and 3.36 indicate that a parameter change AAj^r 

r%1 /-k -P \ ^ • — 



■ H/L.J 'AJo' 

a 

,8) primarily affect ;*:he overshoot of X5(t). 

X 

The plots of j=1,2,...,8 in 

Fig. 3.37 and 3.38 indicate that a parameter change 
(j=1 , 2 ,. . . , 8) primarily affect the overshoot of X6 (t) . 

The plots of Ayj’ = j=1»2,...,8 in 

Fig. 3.39 and 3.40 indicate that i parameter change 
(j= 1 » 2 ,. . . , 8) primarily affect the overshoot of X7(t). 

The plots of /Afij’ = j=1,2,...,8 in 

Pig. 3.41 and 3.42 indicate that a parameter change 
(j = 1 » 2 ,. . . , 8) primarily affect the overshoot of X8(t). 

The plots of AgJ = j=1,2,...,8 in 

Fig. 3.43 and 3.44 indicate that a parameter change A AJq, 

(j=1 » 2 ,. .. , 8) primarily affect the overshoot and rise time 

of X9(t) . ^ ^ 

V U X j o 

The plots of A|oj = j=1,2,...,8 in 

Fig. 3.45 and 3.46 indicate that a parameter change 
(j=1 » 2,. .. , 8) primarily affect the rise time of Xl0(t). 
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The plots of 



^ X II 

‘‘U ■ j-lr2,...,8 in 

Fig. 3.47 and 3.48 indicate that a parameter change 2iA3o, 
( 3 =“* » 2 ,. . . , 8) primarily affect^he overshoot of lC11(t). 

The plots of /\ ( 2 ^ “ “ I /\ j 3“^»2»...r8 in 

Fig. 3.49 and 3.50 indicate that a parameter change 
(j=1 , 2 , 8) primarily affect the overshoot of X12(t). 

Table II and III show the above analysis in a 
condensed way to give a general picture of all states and 
output sensitivity functions with the correspondent effect 
as function of time. 

2. Parame t er- In duced O utpu t An al ysi s 

As Shown previously, if ^ oc << , the actual 

output can be written as 



y(t,oc) = y(t , 



+ ^ (t, o^) ^ 



(3.34) 



The computer program given in appendix E was written 
for simulating the system when 10??, 30%, and 43% of varia- 

tion from the nominal value of each parameter is 
assumed. Fig. 3. 51 through 3.62 give the actual output when 
each parameter is varied 10% from the nominal value. 

Fig. 3. 63 and 3.64 give the actual output for 30% of 
variation from the nominal value of each parameter assumed. 
Fig. 3.65 and 3.66 give the actual output for 43% of varia- 
tion from the nominal value of each parameter assumed. 
These plots show the output of the state variables X3 and 
XII that present strong deviations just to give the behavior 
of the system when parameter variations are not small. One 
notes that modeling is starting to break down. 
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Figure 3,2 Sensitivity of XI with Hespect to A1,&2^A3,A4. 
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Pigure 3.3 Sensitivity of 12 with Respect to A1,A2^A3^A4. 
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Figure 3,4 Sensitivity of X3 with Respect to &1,A2,A3,A4. 
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Figure 3.5 Sensitivity of X4 with Respect to A1,A2,A3,A4. 
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Figure 3.6 Sensitivity of X5 with Respect to A1,A2,A3,A4. 
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Figure 3.8 &ctuei and Noeinal Output of XI (10X variation) • 
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Figure 3,9 



ictual and Nominal Output of 12 (10% variation) 
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Figure 3.10 



Actual and Nominal Output of X3 



(10% variation) . 
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Figure 3.11 



Actual and Noiinal Output of X4 



(10% variation) . 
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Figure 3.12 Actual and Nominal Output of X5 (10?5 variation) . 
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Figure 3.13 Actual and Nominal Output of X6 (10% variation) . 
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Figure 3.14 Actual and Nominal Output of XI (33% variation). 
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Figure 3.15 Actual and Nominal output of X2 (33% variation) . 
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Figure 3.16 



Actual and Nominal Output of X3 



(30 X variation) . 
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Figure 3.17 Actual and Nominal Output of 14 (30 variation) . 
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Figure 3.18 



Actual and Noainal Output of 15 



(30% variation) . 
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Figure 3.19 



Actual and Nominal Output of Z6 



(33^ variation) . 
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Figure 3*20 Actual and Noainal Output of XI (43% variation) 
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Figure 3.21 



actual and Nominal Output of X2 



(49% variation) . 
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Figure 3.22 



Actual and Nominal Output of X3 



(43 36 variation) . 



74 



i 



I 

I 



on 



Ok 



Figure 3.. 




23 Actual and Noninal Output of X4 (40% variation) . 
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Figure 3.24 



Actual and Nominal Output o£ X5 



(49 X variation) . 
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Figure 3.25 



Actual and Nooinal Output of X6 (43% variation) . 
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Figure 3.26 Boll-Yav Nominal and Sensitivity Hodels. 
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3.27 Sensitivity of XI with Respect to h1,A2,A3, 44 
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Figure 3.28 Sensitivity of XI with Respect to X5,A6,A7, A8. 
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Figure 3.29 Sensitivity of X2 with Respect to MrA2,A3,a4. 
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Figure 3.30 Sensitivity of X2 with Respect to &5,46,A7,A8. 
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Figure 3.31 Sensitivity of X3 with Respect to &1,A2,A3,A4. 
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Figure 3.32 Sensitivity of 13 witli Respect to &5,A6,A7, A8. 
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Figure 3.33 Sensitivity of X4 with Respect to M,A2,A3,A4. 
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Figure 3.34 Sensitivity of X4 with Respect to \5,A6,A7 ,a8. 
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Figure 3.35 Sensitivity of X5 with Respect to M,A2,A3,A4. 
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Figure 3.36 Sensitiyity of X5 with Respect to h5,A6,A7, A8. 
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Figure 3,37 Sensitivity of X6 with Respect to hi ,12,43, A4. 
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Figure 3.38 Sensitivity of 16 with Respect to 45,A6,A7,A8. 
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Figure 3.39 Sensitivity of X7 with Bespect to &1«A2,A3,A4. 
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Figure 3.40 Sensitivity of 17 witli Respect to A6, A7, A8. 
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Figure 3,41 sensitivity of X8 with Respect to ai,a2,A3eA4. 
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Figure 3*42 Sensitivity of X8 witli Respect to &5,A6,A7,A8, 



94 




I 



\ 



i 

i 

i 

( 



I 




Figure 3.43 Sensitivity of X9 vitti Sespect to &1,A2,A3,A4. 
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Figure 3.44 Sensitivity of X9 with Respect to &5,A6,A7,A8. 
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Figure 3.45 Sensitivity of X 10 with Respect to i 1 ,A 2 ,A 3 ,A 4 . 
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Figure 3.47 Sensitivity of XII with Respect to A1,A2,A3,A4 



99 



1 






I 

I 



i 



LEGEND 







3aa 

00-s7 QQ*ds 00*S2 00^^ QQ-SZ- QQ-QS- QO St- 00 -001- 00 S2I- 

sin 



00*sl 00*09 OQ'Sh 00“0E 00* SI 00*0 00 ’SI- 00*^- 00*Sfl- 

9in 



00-Sl 05*21 00*01 



OS *4 00 *S OS *2 

LUl 



00*0 QS*2- Q0*s- 



00*09 00*0b 00*02 






00*0 00*02- 00*0h- 00*09- 00*08- 00*001- 

9in 



I 

I 



Figure 3.48 Sensitivity of X11 with Respect to A5,A6,A7,A8, 



100 



LEGEND 




pen 

QS*0 q5tQ oPq 02*0 oro OO^Q oro- 02’d- oFo- 



9f0 2TQ 80*0 



h0*0 00*0 fi0*0- 80*0- 

2211 



2fd- 91 *0- 



26*0 h2*0 Sp^ 



80*0 00*0 80*0- 91*0- 

£211 



1i2*Q- 2£’0- 



9t*o iro 80*0 ho*o 00*0 iio*d- 80*d- 2i*d- 9i*d- 

fi2I1 



Figure 3.49 Sensitivity of X12 with Respect to A1rA2rA3^A4, 



101 



LEGEND 







Pigure 3.50 Sensitivity of X12 with Respect to A5,A6,A7,A8. 



102 







TABLE II 

Influence of Parameters 






UJ 

2 




m 

Z 




1 


Ul 

U 


US 

-1 


ui 

Z 






Ul 

J 


a 


% 




N 


UJ 

Z 


Ul 


lii 

r 






UJ 

Z 


01 

J 


Ul 

z 






Ui 


UJ 


Ui 

2 




Vo 


Ui 

Z 


UJ 


Ul 

I 


Vt> 

io 

x: 


UJ 

Z 


Ul 

u 


Ul 

£ 




vD 


Ui 


UJ 

-J 


Ui 

2 




uy 


UJ 

Z 


a 


ul 

X 


IC 

tp 

/C 


uJ 

Z 


UJ 

-I 


Ui 

X 




in 


UJ 

J 


UJ 

J 


10 

Z 


' 


▼* 


UJ 

z 


UJ 


lu 

1 


5 

/<: 


111 

2 


Ul 

J 


UJ 

2 




UJ 

2 


Ul 

-J 


LiJ 

X 


«o 


01 

u 


UJ 

-1 


Ul 

z 


»o 

v> 

/< 


Ul 

-i 


Ul 


UJ 

X 


lO 

KO 

/< 


Ul 

J 


UJ 

J 


UJ 

z , 




rvi 


UJ 

Z 


UJ 

-1 


UJ 

z 


rJ 

to 

/< 


Ul 

r 


Ul 


UJ 

2 




r>* 

sO 

/< 


UJ 

2 


UJ 


Uj 

X 


1 

1 




UJ 

r 


Ul 

J 


!iJ 

z 


'5 

/< 


Ul 

z 


UJ 

_J 


UJ 

2 




J 

/< 


UJ 


Ui 

-J 


ui 

2 




Ut ui 
^ K 


V- 

o 

0 

1 
CO 

a: 

Ul 

> 

a 


y 

o ui 

Ul < 

y h 




UlUl 

dh 


H 

o 

o 

X 

lO 

Di 

Ui 

> 

o 


y 

0 id 

4 y 

"5 

iO c/1 




|UI 

Sz 

5y 


y 

o 

0 
X 

U9 

0^ 

01 
> 
0 


^ 01 
UJ ^ 

y 

Sen 




















C 


UI 

J 


ul 

iD 


UJ 

Z. 






!3 


m 

cJ> 


UJ 




00 

jC 


in 

J 


Ul 

tO 


UJ 

X 


■ 

I 


N. 

/< 


UI 

«j 


ul 

U) 


ui 

z 


K 

CM 

x: 


Ul 

2 


Ul 

i/i 


UJ 

X 




Ul 

2 


lu 

w 


Ui 

2 ; 


<< 


UJ 

u 


ul 

U> 


Ul 

z 


«M 


UJ 

-J 


UJ 


U) 

X 


VP 

»n 

/< 


Ul 

z 


Ul 

0 


UJ 

Z 


lO 

/< 


UJ 


Ul 

\0 


Ul 

z 


to 

CM 

/C 


UJ 

«J 


Ul 


Ul 

X 


If* 

ro 

K 


ai 

Z 


UJ 

w 


U) 

Z 




r 

y< 


UJ 


ui 

(0 


Ul 

X 




/<: 


Ul 

X 


Ul 

SP 


UJ 

X 


K 


Ul 

2 


UJ 

V3 


Ul 

2 




51 

/< 


UI 


ul 

CO 


lU 

2 




fo 

CM 

/< 


UJ 

-1 


Ul 

lO 


UJ 

2 




in 

/C 


Ul 


UJ 

to 


Ul 

2 




r** 

/C 


ul 

z 


ut 


Ul 

2 


«si 

d 

/< 


Ul 

2 


Ul 

\p 


Its 

X 




fM 

»n 

/< 


Ul 

2 


Ul 

to 


Ul 

Z 




/< 


UJ 

r 


UJ 

iO 


Ui 

2 




N 

/< 


UJ 

2 


Ul 

o 


iJ 

z 




!n 

/< 


Ul 

X 


Ul 

to 


Ul 

Z 






lllUl 

2 H 


V- 

o 

0 

1 
CO 

Ul 

> 

o 


y 

O Ui 

^1: 
Ul < 

V- y 

vO^ 




UJU^ 
jo 2 

5y 


y 

o 

o 

X 

J) 

Di 

UJ 

> 

o 


y 

0 id 

4 y 

y *: 

cO 






5y 


y 

o 

0 
X 
to 

01 
> 
Q 


P 01 

<b 

lij < 

H y 
Sen 






X 













103 



I 



I 





Influence of 


TABLE III 

Parameters (cont. Table III 




<k> 

o 

/C 


m 

j 


Ul 

-J 


Ul 

2 




<o 

/< 


Ul 

-1 


Ui 


ui 

Z 




CD 


iiJ 

X 


Ui 

t/7 


Ui 

z 




o 


ui 

•/ 


ul 


01 

z 


/< 


01 

2 


Ui 


Oi 

z 




rv 

/< 


Ul 

2 


Ul 

lO 


ul 

2 




>s 

o 


!3 


in 

-1 


ul 

2 


NJ 

/<* 


Ul 

2 


Ui 


z 


VP 

/< 


Ul 

-J 


Ul 

(/) 


Mi 

2 




O 

/< 


UI 

U 


(il 

-1 


m 

2 


V» 

/C 


Ul 

z 


ui 


(il 

z 




M 

K" 


Ui 

z 


Ui 

V) 


Ul 

2 




•V 

o 


Ui 


UJ 


Ul 

2 


/< 


ul 

2 


Ul 

cl) 


Jii 

z 




«<■ 


Ul 

2 


Ui 

w 


Ui 

2 


o 


a 


a 


ill 

2 


— 

/<' 


Ui 

J 


Ml 

V) 


<u 

z 




/<“ 


Ui 

J 


Ui 

o 


Ui 

2 


o 

/< 


UI 


HI 


UJ 

z 


INI 

/< 


lij 

2 


ul 

tl> 


0/ 

z 




»%* 

/T 


Ui 

2 


Ui 

U7 


Ui 

2 




o 

/< 




UJ 


UJ 

2 




Ul 

J 


Ul 


Ul 

2 


<»* 


'Oi 

2 


Ui 

Ul 


Ui 

z 




a|U> 
Ui Z 

5h 


V- 

o 

o 

X 

(/) 

a: 

Ul 

>. 

a 


y 

o ul 

Ui ^ 

UH 

\n <n 






UlUl 

a? 


H 

O 

0 
2 

01 
Ci! 
UJ 
> 
a 


> 

0 uJ 
A h 

</l *il 




iUl 

Sz 

cdh 


OVERSHOOT 


STEADY 

STATE 






o 

X” 




x" 






CM 

X“ 


' 




00 

/C 


ul 

J 


a 


itf 

I 




<0 


ui 

2 


ui 

1/9 


ul 

2 




cO 

« 


a 


Ul 

-J 


Ui 

z 


»». 

K 


UI 

Z 


a 


(J 

2 


CD 

/< 


Ui 

Z 


Ul 


Ul 

z 




t- 

CP 


Ul 

-1 


a 


lU 

i 


ii 

/< 


UI 

J 


ui 


111 

2 


\s 

cO 

x: 




Ui 

0> 


01 

2 


(P 

/< 




Ui 

-1 


Ui 

Z 


/< 


01 

2 


lU 

J 


Ul 

Z 


00 

/<: 


ui 

2 


Ui 

lO 


Ul 

2 


w-t 


u; 

-J 


Ul 

J 


Ul 

2 


? 

/< 


111 

2 


Ul 

J 


Ul 

X 


'T 


ni 

Z 


UJ 

4> 


Oi 

z 




Ui 

-J 


Ui 

-i 


ui 

X 


lO 

e- 

/< 


UI 

J 


, a 


UJ 

2 


•O 

<» 

/< 


til 

2 


Ul 


Ui 

z 




Ul 


Ui 

-J 


Ui 

2 


W 

/< 


UI 

2 


a 


01 

2 


M 

01 


ni 

:t 


Ul 

lO 


al 

2 




ff 

/< 


Ui 

-1 


Ui 


in 

2 


/< 


UJ 

-J 


111 


ul 

2 


£ 

K 


Ul 

2 


Ul 


Ui 

z 




s 

/< 


Ul 


In 

_i 


Ul 

z 






(ItUl 
ui Z 


w 

o 

o 

X 

(0 

(t: 

Ui 

> 

a 


y 

oui 

UJ < 

W- H 

\n <n 






ujui 
12 ^ 


H 

o 

0 
2 

01 
Di 

Ui 

> 

a 


y 

0 si 
4 y 
«j< 

y 






lid 

a? 

5h 


1- 

o 

0 

1 

U9 

U 

Ui 

> 

Q 


O (11 

<b 

UJ < 

5icn 








OD 

X 






x*^ 


> 



104 



r 



00 - 



00 - 



I 

Figure 




oas/3aa 



le 00'9t OO'Sl 00*81 



00*6 00*9 

SIX 



00 •£ 00*0 00*E- 



3.51 Actual and Nominal Output of XI (10J6 variation) 
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Figure 3.52 Actual and Nominal Output of X2(1035 variation). 
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Pigure 3.53 Actual and Nominal Output of X3(10* variation). 
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Figure 3.54 actual and Nominal Output of X4 (10X variation) . 
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Figure 3.55 Actual and Nominal Output of X5(10A yariation) 
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Pigure 3.57 Actual and Nominal Output of X7(10X variation). 
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Figure 3.58 Actual and Nominal Output of X8(10X variation). 
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Figure 3.59 Actual and Nominal Output of X9(10X variation). 
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Figure 3.60 actual and Nominal Output of X10(13* variation). 
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Figure 3.62 Actual and Nominal Output of X12(1()% variation). 
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Figure 3.63 Actual and Nominal Output of 13 (30)6 variation) . 
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Figure 3.64 Actual and Noninal Output of XII (33 X variation). 



118 




119 



TIME 




Figure 3.66 Actual and Nominal Output of 111(43% variation). 
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17. ftP PLICATIO H OF S ENSITI VITY ANAtXSIS TO NO^INEAR SYSTEMS 

A. INTRODDCTION 

The three dimensional nonlinear system used here for 
sensitivity analysis purposes is presented in appendix c. 
The missile configuration, size and mass properties are the 
same as those used in the lineat system and were presented 
in appendix A, 

In section B the sensitivity equations with respect to 
parameters of nominal system is shown and in section C the 
sensitivity equations with respect to parameters that showed 
to be of more effect in the previous analysis is presented. 

Section D shows the results of the trajectory sensi- 
tivity equations when step inputs are applied at specific 
time as mentioned in appendix C. 

Section E gives the parameter-induced output analysis 
when each parameter is varied 10% from the nominal value. 

B. NONLINEAR EQUATIONS OF THE NOMINAL SYSTEM 

From the block diagram of Fig. C.1 and Egn. C.3 through 

C. 21 as given in appendix C one have the following nominal 
equations. 



XI = - C4 XI - C02 C4 A2 (4.1) 



X2 = - C5 X2 - C6 NZC ♦ C6 C19 Cos (X7) + C7 XI (4.2) 
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X3 = (C5 C8 - C9) X2 + C6 C8 NZC (4.3) 

- C6 C8 C19 Cos(X7) - C7 C8 XI + (C8/Conv2) XI 7 X18 
+ Cl C8 A1 + (C9/Conv) X5 

X4 = - C3 X4 + C3 Conv X3 (4.4) 

X5 = (X17 X18)/Conv + Cl Conv A1 (4.5) 

X6 = X5 - KB C02 A2 - (X16 X18)/Conv (4.6) 

X7 = X5 Cos(X19) /Conv - X17 Sin (XI 9) /Conv (4.7) 

X8 = - CIO X8 ♦ (C12 - Cl 1 C13) XIO (4.8) 

+ C11 C14 PHC - C11 C14 X19 - C11 CO (A8 X16 + A7 X15 
♦ A3 X13) - (C12/Conv) X18 

X9 = - X9/T1 + K1/T1 C02(AA X16 ♦ A6 X15) (4.9) 

4 

XIO = - C13 XIO + C14 PHC - C14 X19 (4.10) 

XII = - C15 XII + C16 C0(A8 X16 ♦ A7 X15 + A3 X 13) (4.11) 



122 



X12 = - C17 X12 + K(C17 - (C10/C18) X8 



(4. 12) 



(K/C18) (C12 - C11 C13) X10 ♦ K Cl 1 (C14/C18) PHC 

- K C11 (C14/C18) X19 - (K/C18) CO A8(C11 - (K/C18) 

CO A3(C11 *■ C16) X13 - (K/C12) (Cl 8 Conv) X18 

+ K(C15/C18) - C17) XII 

X13 = - C3 X13 ♦ C3 Conv X12 (4.13) 



X14 = (K2/10) (- (X9/T1) C02(AA X16 ♦ A6 X15) (4.14) 

- X15 (X18/Conv2) + C01(A9 X16 ♦ A4 X15 ♦ A5 X13) 

- (KYP/Conv2) (X5 - KB CO 2 A2 - XI 6 (X 1 8/Con v) X18 

- (KYP/Conv) XX6 C0(a8 X16 ♦ A7 X15 + A3 X13) 

♦ (K2/Conv) X17 - K2 (KYP/Conv2) XX6 X18 

X15 = - C3 X15 ♦ C3 Conv X14 (^.15) 



X16 = KB C02(AA X16 ♦ A6 X15) ♦ X5(X18/Conv) - X17(4.16) 

X17 = - (X5 X18)/Conv + C01 Conv(A9 X16 (4.17) 

+ A4 X15 + A5 X13) 

X18 = CO Conv(A8 X16 + A7 X15 A3 X13) (4.18) 



XI 9 = Xl8/Conv 



(4. 19) 



The correspondence of nonlinear state vectocsi as given 
in the above equations with relation to the nonlinear system 
presented in Fig.C. 1 in appendix c are : 

XI = X , X2 = Y , X3 = , X4 = Sp , X5 = q , 

X6 = , X7 = © , X8 = Y1 , X9 = Y2 , XI 0 = XI , 

XII = X2 , X12 = Src» X13 = Sr, X14 = X15 = 8^, 

X16 = , X17 = r , X18 = p , and X19 = 4 • 



The parameters of interest for this nonlinear system are 
given by 

A 1 a2 =C}<X,ip\, A3 =C„U), a 4 = , and A5 



M 



"«Y 

Definition of the constants CO through C19 are given in 
appendix C. 

As previously, this nomenclature is used here to easily 
apply the sensitivity theory. 



C. HOHLIHBAR SENSITIVITY EQUATIONS 

From the sensitivity theory given in chapter 2 one knows 
from Egn.2.12 that for nonlinear systems one has 

'1)X 

^ = — + » = 0 (4.20) 

Here 



X = f {X1,X2,.. .,X19,t»NZC,Al,A2,... ,A5} 



(4.21) 



iThe state vectors and parameters in the nonlinear case 
have no correspondence to the stats vectors and parameters 
used in the linear case. The correspondence can be found in 
the appendices B and C. 
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From Eqn.2.14 one knows that the trajectory sensitivity 
vector is given as 





-li 



C<0 



A (0) = 0 



(4.22) 



Applying the above theory one obtains the nonlinear 
sensitivity equations as given in Eqn.4.24 through 54e42, 
where ^ was considered the parameter to be varied. 



Ill = - C4/111 - C02 C4{DA26^61 ♦ DA24a41) 



(4. 23) 



A21 = - C5/121 - C6 C19 Sin(X7) A71 + C7^11 (4.24) 



^31 = (C5 C8 - C9) J.21 + C6 C8 C19 Sin (X7) ^71 (4.25) 

- C7 C8 All + C8(x18A171 + (X17 A 18 1 ) / (ConvZ) 

+ Cl C8(DA16A61 +DA14A41 + 1.) + (C9A51)/Conv 

A41 = - C3A41 + C3 Conv A3 1 (4.26) 



A51 = (X17A181 ♦ X18A17 1)/Conv (4.27) 

+ Cl Conv(DA16A61 + DA14A41 + 1.) 

A61 =A51 - KB C02(DA26A61 + DA24A41) (4.28) 

- ((X16A181 + X18 A161) /Conv 
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.171 = (Cos(X19) /i51 - X5 Sin (X19) .1 191) /COQV (4.29) 

- (Sin(Xl9) ^171 + X17 Cos (X19) ai71) /conv 

Jlai = - CIO .181 + (C12 - C11 C13) A^0^ (4.30) 

- C11 C14.1191 - C11 CO (A8/1161 ♦ A7.1151 + 43^131 

+ /161 (XI 6 Da8 ♦ X15 Da7 +X13 DA13)) - (Cl 2 .11 81) /Conv 

yL9^ = -.191/T1 ♦ K1/(T1 C02) (AAxl161 (4.31) 

+ A6 /1151 *■ A6\ (X16 DAA + X15 Da6) ) 

# 

/1101 = - C13.1101 - C14/1191 (4.32) 

/ 

m 

.1111 = - C15/1111 + C16 C0(A8.1161 + A7^151 (4.33) 

+ A3/1131 ♦ .161 (X16 DA8 ♦ X15 DA7 + X13 DA3) ) 

A121 = - C17.^121 + K(C17 - C10/C18) .181 (4.34) 

+ (K/C18) (C12 - C11 C13) .1101 - K C11 (C14/C18) ^191 

- (K/C18) C0(C11 ♦ C16) (A8.1161 + A7.1151 + A3 .1131 
♦ A3 4.131 + 4.61 (X16 DA8 ♦ X15 DA7 + X13 Da3) ) 

- K C12/(C18 Conv) 4.181 ♦ K (C15/C18 - C17)4111 

4.131 = - C3 4.131 + C3 Conv4121 (4.35) 



126 



^141 = (K2/10) (- X9/T1 + (K1/T1) C02(AA 



(4.36) 



/i161 + A6/1151 + 261(X16 DAA + X15 DA6) ) - (X15/1181 
+ X18 Z151) / (Conv2) + C01(A9.1161 + A4^151 + A5 
+ ^61(X16DA9 + X15 DA4 + X13 DA5) ) - KYP/(Coav2) 

(X5 2181 ♦ X18/151 - KB C02(A2^181 + X18(DA25^61 
♦ DA24/141) -{X182ai61 + 2 X18 X16/1181) /Conv)) 

- (KYP/Con\r) C0{ xx6/i161 + X16ai61) + A7ai61) 

+ A7(XX6/1151 ♦ X15/1161) + A3(XX6/i131 + X13/1161) 

+ X16(XX6 DA8/161 + A8^161) + X15{XX6 DA7xl61 
+ A7 /1161) + X13(XX6 DA3 ^61 + A3 ^161))) 

+ (K2/Conv) /1 171 + K2/191 - K2 KYP/(Conv2) 

(XX6 /1181 ♦ X18 ^161) ) 

^151 = - C3/1151 + C3 Coavi.141 (4.37) 



A161 = KB C02(AA ,161 + A6 2151 +261( (4.38) 

X16 DAA + X15 DA6)) + (X62181 + X 1 8 26 1) /Conv -2171 

2171 = - (X5 2181 + X18251)/Conv (4.39) 

+ C01 Conv(A92161 + A4 2151 + A5 2131 +261(X16 DA9 
+ X15 DA4 ♦ X13 DAS) ) 



2181 = CO conv (A8 2161 + A7 2151 



(4.40) 



+ A3 2131 + 261(X16 DA8 ♦ X15 DA7 + X13 DA3)) 



2191 =2181/Conv 



(4. 41) 
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The same procedure above must be made for the other 
parameters that most affected the previous analysis and will 
be given in the nonlinear sensitivity analysis presented in 
the next section. 



the nominal 
Here, the 
cases given 
in the time 



D. NONLIHEAR SENSITIVITY ANALYSIS 

A computer program in appendix 3 shows the simulation of 
the sensitivity eguations with respect to the selected 
parameters. 

The number of equations solved are 19 for 
equations and 95 for the sensitivity equations, 
parameters were selected from the previous two 
in Chapter III that showed to have most effect 
response of the systems. 

From the analysis performed, two paraueters were 
selected from the uncoupled pitch autopilot and three from 
the coupled roll-yaw autopilot, which has been shown to be 
more sensible to parameter variations. These parameters are 
respectively (A2) and Cff^^ (A4) in the uncoupled pitch 

autopilot and ( A3) , (A5) , and in the 

coupled roll-yaw autopilot. 

In order to compare the influence of the parameter of 
interest in the nonlinear case some state variables were 
chosen to be analysed. The state variables selected were: 

Prom the uncoupled pitch autopilot : 

X3 = Spg , X4 = Jp , X5 = q , X6 = X 

From the coupled roll-yaw autopilot : 



X12 = <S 



Rc 



X16 = 



, X13 = , X14 = , X15 = Sy 

, X17 = r , X18 = p , and X19 = ^ . 
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The results are plotted in Fig. 4.1 through t 4.12 and 
Table IV and V. Each state variable output is plotted 
separately with the correspondent five output sensitivity 
functions. 

From the plots, the following observations can be made; 

Fig. 4.1 indicates that the state variable X3(5f>^) is 

little affected in the rise time by the parameters A1 and 
A2. The overshoot is strongly affected by A5, A1 and A4 and 
little affected by A2 and A3. The steady state is strongly 
affected by A1 and little affected by A2. 

Fig. 4. 2 indicates that the state variable X4(<Sp) is 

little affected in the rise time by A1 and A2. The over- 
shoot is strongly affected by A1, A3, and A4 with little 

effect by A2 and A5 . The steady state is strongly affected 
by A1 and little affected by A2. 

Fig. 4. 3 indicates that the state variable X5 ( cj ) is 

strongly affected in the rise time by A1 and A2. The over- 
shoot is strongly affected by A3, A4, A5, and little 

affected by A1 and A2. The steady state is not affected by 

parameter variations. 

Fig. 4. 4 indicates that the state variable X6 ( ) is 

strongly affected in the rise time by A1 and A2. The over- 
shoot is strongly affected by A3, A4, and A5. The steady 

state is little affected by A1 and A2. 

Fig. 4. 5 indicates that the rise time and steady state of 
the state variable X12(5f^^J are not affected by parameter 
variations. The overshoot is strongly affected by A3 and A5 
and little affected by A1, A2, and A4. 

Fig. 4.6 indicates that the rise time and steady state of 
the state variable X13(Jrj^) are not affected by parameter 
variations. The overshoot is strongly affected by A3, A4, 

and A5 and little affected by A1 and A2. 

Fig. 4. 7 indicates that the rise time and the steady 
state of the state variable X14 (Sy^) are not affected by the 
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parameter variations. The overshoot is strongly affected by 
A4 and A5 and little affected by A1, A2, and A3. 

Fig. 4. 8 indicates that the rise time and the steady 

state of the state variable X15(6y > are not affected by the 
parameter variations. The overshoot is strongly affected by 
A3, A4, and A5 and little affected by Al, and A2. 

Fig. 4. 9 indicates that the rise time and the steady 

state of the state variable X16(^ ) are not affected by the 
parameter variations. The overshoot is strongly affected by 
A3, A4, and A5 and little affected by Al, and A2. 

Fig. 4. 10 indicates that the rise time and the steady 

state of the state variable X17( ) are not affected by the 

parameter variations. The overshoot is strongly affected by 
A3, A4, and A5 and little affected by Al, and A2. 

Fig. 4. 11 indicates that the rise time and the steady 

state of the state variable X18( ( are not affected by the 

parameter variations. The overshoot is strongly affected by 
A3, A4, and A5 and little affected by Al, and A2. 

Fig. 4. 12 indicates that the rise time and the steady 

state of the state variable X19( (j> ) are not affected by the 
parameter variations. The overshoot is strongly affected by 
A3, A4, and A5 and little affected by Al, and A2. 

E. PAHAHETER-INDDCED OOTPOT ANALYSIS 

As shown previousy, if Aoi« , the actual output can 
be written as 



y (t , «< ) = y (t, ) ♦ (T (t , (y^ ) 2 :^ oc 



(4.42) 



The computer program given in appendix G was written for 
simulating the system when variations from the nominal value 
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of each parameter is assumed. Figs. 4, 13 through t 4.24 give 
the plots of the actual and nominal output of the state 
variables that are present in the unclouped pitch autopilot 
and coupled roll-yaw autopilot when 10S of parameter varia- 
tion is assumed. The others state variables are not given 
here because they are just intermediate state variables. 

these results confirm the theory presented in chapter 2 
when small parameter variations are assumed. Figs. 4.25 and 
4.26 present the actual and nominal output of the state 
variables X3 and XII that showed pronounced variations when 
30K of parameter variation is assumed. Figs. 4.27 and 4.28 
present the actual and nominal output of X3 and XII when 40% 
of parameter variation is assumed. From these polts one 
note that the modeling is starting to break down. 
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Figure 4.1 Sensitivity of X3 with respect to &1 , A2«A3,A4, AS. 
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Figure 4.2 Sensitivity of X4 with respect to A3 ,A4,A5 . 
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Figure 4.3 Sensitivity of X5 with respect to &1 ,&2,A3,A4,A5 
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Figure 4.4 Sensitivity of X6 with respect to 41 ,A2,A3,A4, A5- 
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Figure 4.5 sensitivity of XI 2 with respect to i 1 , A2, &3, A4, A5. 
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Sensitivity of XI 3 with respect to M,A2, A3, A4,i5. 
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Figure 4.7 Sensitivity of X14 with respect to M,A2, A3,A4, &5. 
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Figure 4*8 Sensitivity af X15 with respect to 1 1 ^ A2^ A3, A4^ AS. 
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Figure 4,9 Sensitivity of X16 with respect to & 1 A3, A4, &5. 
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Sensitivity of X17 with respect to A1,A2, A3,A4 , a5. 
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Figure 4.11 
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Sensitivity of X18 with respect to A1» A2, A3,A4^ A5. 
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Sensitivity of X19 with respect to Al» A2,A3»i4, A5. 
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Figure 4.14 Actual and Nominal Output of X4 (10X variation) . 
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Figure 4.15 actual and Hoainal Output of X5 (10 S variation) . 
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Figure 4.16 Actual and Noainal Output of X6 (13% variation). 
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Figure 4.17 Actual and Noninal Output of XI 2 (10X variation). 
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Figure 4*18 Actual and Nominal Output of XI 3 (10X variation) 



149 






Figure 



□ G 

Q 

Z 

UJ 



uj cn 





II 



2f^ BO*Q 90*0 eO‘Q '□□‘ll EO'Q- 30‘Q- 60‘d- 2T0- 

ftU 

pen 



aro BQ*0 90*0 E0*0 00*0 E0*0- 90*0- 60*0- 21*0- 

SliU 



- J 

4,19 Actaal and Noainal Output of X14 (10% variation)* 



150 



'oVoO iVoO 2.00 3.00 «i'.00 s'. 00 oToO 7.00 Q. 00 fiLOO To . 00 

TIME 

Sec 







":ll 



1 



□ o 



a 

z 

LU 




O 

O 




00 ■3 TE'K E9-'z fipo SfO- Hh"^- El'ji- tS*S- 0S*i- 

SSU 



Figure 4.20 Actual and Nominal Output of X15 (10% variation). 
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Figure 4.21 Actual and Honinal Output of X16 (10* variation). 
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Figure 4*22 Actual and Nominal Output of 117 (10% variation). 
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.23 Actual and Noainal Output of X18 (10% variation). 
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4.24 Actual and Nominal Output of XI 9 (10% variation). 
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Figure 4.25 Actual and Noninal Output of X6 (30 X variation) . 
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4^27 Actual and Nominal Output of X6 (43% variation) 
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Figure 4.28 Actual and Honinal Output of 116 (40% variation). 
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V, CONC LDSION S AND RBC3MME NDAT IOHS 



A. COHCLUSIOHS 

The results of the nonlinear 3-D parameter sensitivity 
analysis presented in chapter IV verify rhe linear analysis 
of chapter III. These conclusioas can be quickly made, 
using Tables I, II, III and IV that give a brief review of 
the time response of the linear and nonlinear system. 
Another means of ccmparison of the linear and nonlinear 
analysis is to use the figures that give a precise visuali- 
zation of the time response. 

The comparison of the linear and nonlinear analysis can 
be done by means of the plots as follows: 



B. EFFECT ON STATE VABIABLES DDE TO PARAMETER VARIATIONS 

1 . Effec t on Sta te Varia ble D ue to Var iati on in 

tn 

Figs. 3. 7 and 4.1 ( Ac 4 and ) show that is 

little affected in the rise time, strongly affected in the 
overshoot and steady state. 

Figs. 3. 4 and 4.2 ( and X 4 j ) show that 6 f> is 

little affected in the rise time and strongly affected in 
the overshoot and steady state. 

Figs. 3. 2 and 4.3 (X,^ andAs”!) show that q is 

strongly affected in the rise time, little affected in the 
overshoot and not affected in the steady stats. 

Figs. 3. 3 and 4.4 ( A^l ^ show that «< is 

strongly affected in the rise time, little affected in the 
overshoot and steady state. 
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2 . 



Effect on Sta te Variable Du e to Va r iati on in 

C_(ot,iE) : 

M 





Figs. 3. 7 


and 


( XfcZ A 32 ) show 


t hat 


^Pc is 


little 


affected 


in 


the rise time, overshoot 


and 


steady 


stat e. 




Figs. 3 


.4 


and 4.2 (^ 92 . ) show 


that 


Sp is 


little 


affected 


in 


the rise time, overshoot 


and 


steady 


state. 




Figs. 3. 2 


and 4.3 (A 12 and show 


that g is 



strongly affected in the rise time and little affected in 
the overshoot and steady state. 

Figs. 3. 3 and 4.4 (^22 >^€*2) show that is 

strongly affected in the rise time and little affected in 
the overshoot and steady state. 

3. Effect on State Variable Due to Variatioi in C, 

Figs. 3. 40 and 4.6 (A 7 g and Ai23) show that is 

little affected in the rise time, strongly affected in the 
overshoot and not affected in the steady state. 

Figs. 3. 42 and 4.7 (Ag^and A 155 ) show that is not 
affected in the rise time and steady state, and strongly 
affected in the overshoot. 

Figs. 3. 46 and 4.8 ( Aio© ) show that 6 yc, is 

little affected in the rise time, little affected in the 
overshoot and not affected in the steady state. 

Figs. 3. 48 and 4.9 ( Ajig and A 153 ) show that by is 

little effected in the rise time, strongly affected in the 
overshoot and not affected in the steady state. 

Figs.3.32 and 4.10 (A35andA|63 ) show that is 

little affected in the rise time, stronglu affected in the 
overshoot and not affected in the steady state. 

Figs. 3. 28 and 4.11 (A|g and Ai73) show that )r is 
little affected in the rise time, strongly affected in the 
overshoot and not affected in the steady state. 
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Figs. 3. 30 and 4.12 ( A 23 Al63 ) show that p is 

little affected in the rise time, strongly affected in the 
overshoot and not affected in the steady state. 

Figs. 3 . 50 and 4.13 (A( 2 gandAi<j 3 ) show that^is 
not affected in the rise time and steady stats and strongly 
affected in the overshoot. 

4. Effec t on St ate Var iable Due ^ Var iati on in C ; 

* " “ ■ 

Figs. 3. 40 and 4.6 and A| 2 <q) show that is 

not affected in the rise time and steady state and little 

affected in the overshoot. 

Figs. 3. 42 and 4.7 (Ag^ and A 134 ) show that is 

not affected in the rise time and steady state and strongly 

affected in the overshoot. 

Figs. 3. 46 and 4.8 (X{qS and A|.^<^ ) show that is 

little affected in the rise time, strongly affacted in the 
overshoot and not affected in the steady stats. 

Figs. 3. 48 and 4.9 (Aji^ and A^ 5 ^) show that Sy is 
not affected in the rise time and steady state and strongly 
affected in the overshoot. 

Figs. 3. 32 and 4.10 ( As^and A 14”^ ) show that p> is 

little affected in the rise time, strongly affected in the 
overshoot and not affected in the steady stats. 

Figs. 3. 28 and 4.11 ( A 15 and Ajy^) show taat r is not 

affected in the rise time and steady srate and strongly 
affected in the overshoot. 

Figs. 3. 30 and 4.12 (AzS and Aie<^ ) show that p is 
little affected in the rise time, strongly affected in the 
overshoot and not affected in the steady state. 

Figs. 3. 50 and 4.13 (A| 2,5 Aig^) show that ^ is 

not affected in the rise time and steady state and strongly 
affected in the overshoot. 
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5. Effect on St ate Var iabla Dus ^ Variati oi in C : 

nsp. 

Figs. 3. 39 and 4.6 (^kt^^and A|2.6‘) show that S(^<::,is 
little affected in the rise time, strongly affected in the 
overshoot and not affected in the steady state. 

Figs. 3. 41 and 4.7 (^pjand A 135 ) show that is 

little affected in the risr time, stronglu affected in the 
overshoot and not affected in the steady state. 

Figs. 3. 45 and 4.8 (Ai©^ ^ show that 6 yc is 

little affected in the rise time, strongly affected in the 
overshoot and not affected in the steady state. 

Figs. 3. 47 and 4.9 (Aji 3 and Aj^g’) show that Sy is 
little affected in the rise time, strongly affected in the 
overshoot and not affected in the steady state. 

Figs. 3. 31 and 4.10 (A 33 ani A| 6 ^ ) show that is 

little affected in the rise time, stongly affected in the 
overshoot and not affected in the steady state. 

Figs. 3. 27 and 4.11 (A13 and A 17 S ) show that V is 

little affected in the rise time, strongly affected in the 

overshoot and not affected in the steady state. 

Figs. 3. 29 and 4.12 Ajgs) show that p is 

little affected in the rise time, strongly affected in the 

overshoot and not affected in the steady state. 

Figs. 3. 49 and 4.13 (A 123 and Al 3 >j) show that ^ is 
little affected in the rise time, strongly affected in the 
overshoot and not affected in the steady state. 

C. EFFECT OF TRAJECTORY SENSITIVITY FONCTIONS 
1 . Dncoup led Pi tch Aut opilot 

The trajectory sensitivity functions ^11, /i12, 

A13, and in the linear case and the correspondent 

trajectory sensitivity functions 151 and 152 ii the nonli- 
near case show strong effect on the rise time and overshoot 
of the state variable g. 
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The trajectory sensitivity functions yi2 3 and ^24 
in the linear case and the correspondent trajectory sensi- 
tivity function 161 in the nonlinear case show strong 

effect on the rise time of the state variable , 

The trajectory sensitivity functions /133 and ^34 
in the linear case and the correspondent trajectory sensi- 
tivity function /141 in the nonlinear case show strong 

effect on the overshoot and steady state of the srate vari- 
able . 

The trajectory sensitivity functions ^53 and /164 
in the linear case and the correspondent sensitivity trajec- 
tory function ^31 in the nonlinear case show strong effect 
on the overshoot and steady state of the state variable • 

2. Couple d Ho 11 -Yaw Au topil ot 

The trajectory sensitivity functions Z13,/i15, and 
/118 in the linear case and the correspondent trajectory 
sensitivity functions /1175, /11 74, and/i173 in the nonlinear 
case show strong effect on the overshoot of the state vari- 
able r. 

The trajectory sensitivity functions ^23 , /i25, and 
2.28 in the linear case and the correspondent trajectory 

sensitivity functions 2- 185, 2.184, and 2.183 in the nonlinear 
case show strong effect on the overshoot of the state vari- 
able p. 

The trajectory sensitivity functions 2.33, 2.35, and 
238 in the linear case and the correspondent trajectory 

sensitivity functions 2165, >1164, and2.163 in the nonlinear 

case show strong effect on the overshoot of the state vari- 
able ^ . 

The trajectory sensitivity functions 2.73, 275, and 
278 in the linear case and the correspondent trajectory 

sensitivity functions 2 1 25, 2.124, and 2123 in tae nonlinear 
case show little effect on the overshoot of the state vari- 
able 
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The trajectory sensitivity functions /183 , /185, and 
^88 in the linear case and the correspondent trajecrtory 
sensitivity functions d.1 35, /1134, andxl133 in tae nonlinear 
case show strong effect on the overshoot of the state vari- 
able 6 ( 2 . . 

The trajectory sensitivity functions /II 03 and/1105 
in the linear case and the correspondent trajectory sensi- 
tivity functions /1145 and /II 4 4 in the nonlinear case show 
strong effect on the overshoot of the state variable 6^c • 
The trajectory sensitivity function /1108 in the linear case 
and the correspondent trajectory sensitivity fj net ion /1 143 
in the nonlinear case show little effect on the overshoot of 
the state variable Syc. • 

The trajectory sensitivity functions /1 113, /1115, 

and /1118 in the linear case and the correspondent trajec- 
tory sensitivity functions /1 155, /1154, and /1153 in the 

nonlinear case show strong effect on the overshoot of the 
state variable 6 y . 

The trajectory sensitivity functions /i123, A125, 

and /1128 in the linear case and the correspondent trajec- 
tory sensitivity functions /1 195, /1 194, and A193 in the 

nonlinear case show strong effect on the overshoot of the 
state variable . 

The trajectory sensitivity functions have the 
following range of values : 

Minimum Maximum 

Linear - 385.25 632.35 

Nonlinear - 642.9 1772.9 



D. GENERAL CONCLUSIONS 

The parameter C^^ strongly affect the overshoot in 
almost all case. 
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The parameter 
case . 


C., little affect the 

M OL 


overshoot in 


all 


The parameter 
almost all case. 


C«, strongly 


affect 


the 


ov ershoo t 


in 


The parameter 
almost all case. 


Sjp strongly 


affect 


the 


ov ershoot 


in 


The parameter 
almost all case. 


Cp, strongly 


affect 


the 


overshoot 


in 



E. RECOHBENDiTIOHS 

It was not possible to run a CSMP computer program using 
all stare variable (19) and parameters (10) in the nonlinear 
case due to the work area available. The above restriction 

occurred just when the actual and nominal output were 
required to be printed. The computational time in this case 
may be reduced using Fortran subroutines imbedded in the 
CSMP program in order to calculate the necessary parameter 
derivatives . 

In the linear cases the computational time can be 
reduced by means of using the method of Sensitivity Points 
that uses one model for all parameters instead of using one 
model for each parameter. 

The present analysis can be repeated for the case of the 
Circularz airframe given in the [Bef. 2 ], and comparisons 
between both airframes can be performed. 

Future study can be made modeling the system in the 
frequency domain where the "root sensitivity" is analysed. 

Influence of parameter variations on miss distance can 
be analysed using an augumented system and a scenario of 
reference. 



pr 



2The current CSIJP programs in th^ appendices are 
ared ro run the Circular case just by inserting the data 
he correspondent airframe. 
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APPENDIX A 

MISSILE SIZING, HASS PROPERTIES AND AERODYNAIIC DATA 

A. INTRODUCTION 

In this appendix, one considers a model that was 
assumed to be 1/6 scale of the actual elliptical missile 
configuration as given in Fig.A. 1. 

Missile configuration was sized to provide realistc mass 
properties needed for this study. 

Since the main purpose of this work was to perform the 
parameters sensitivity analysis applied to a Sank-ro-Turn 
missile, no effort was expended on a detailed design and 
analysis of -the various parts envolved, which is given in 
[Ref. 2]. Some of figures given in this appendiic were taken 
from the [Ref. 2] fcr easy visualization of the system in 
study. 

B. GEOMETRIC AND MASS PROPERTIES OP MISSILE CONFIGURATION 

One can find in Table VI the properties used in the 
development of the equations for applying the Parameters 
Sensitivity Analysis. 

Fig. A. 2 shows the aerodynamic sign convention, nomencla- 
ture and body-fixed axes where the following three assump- 
tions were made: 

- Plane x^ - zl of Fig.A. 2 is the maneuver plane. 

- Missile is trimmed in pitch (i. e. , My=0 , at fixed 

values of q,c<, and ) . 

-Missile rcll rate is constant. 

Fig. A. 3 shows a block diagram of a BTT autopilot. 
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Figure a.1 Hodel of elliptical coaf iguration ( 1/6 scale ). 
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TABLE VI 

Geonetric and Hass Properties of Missile 
Conf iguratioa. 



WEIGHT (lbs) 24 75 

f't) ^ 

Iyy(slug sq ft) 790 

I Ijj^Cslug sg ft) 3 53 

LENGTH (in) , 1 168 

CENTER OF GRAVITY 100.8(0.6 1' 

distance from nose (in) 

HAX. MAJOR AXIS (in) 41.7 

MAX. MINOR AXIS (in) 13.86 



Inertial acceleration commands are applied in polar coordi- 
nates (i.e. , magnitude of command applied to the 

pitch autopilot and the direction ( 4>c ) applied to the 

roll autopilot) . 

The yaw autopilot is slaved to the roll autopilot to mini- 
mize sideslip angle by coordinating the missile yaw and roll 
motion. 

Achieved maneuver plane or inertial acceleration in rectan- 
gular coordinates (i.e,, and is determined by 

resolving achieved body-fixed accelerations 

(i.e.,"!]^^ and ) through missile roll angle(4> ) (i. e. , 

Euler angles ©• and are assumed to be sufficiently 

small) . 
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Figure A. 2 Aerodynamic Sign Convention. 
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Figure A. 3 BTT AaiOPILOT. 
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APP ENDI X B 
LINEAR SYSTEM MODELS 



A. INTHODOCTION 

The following appendix addresses the informations of the 
linear uncoupled pitch channel autopilot and the roll-yaw 
coupled channel autopilot of a elliptical airframe taken 
from [Ref. 2], 

The linear time domain analysis of the CBTT autopilot was 
assumed that the missile is initially in the desired 
maneuver plane and trimmed at ten degrees angles-of-attack 
(i.e., the equilibrium or trim angle-of-attack in the 

model of Fig.A.I equals 10 degrees and the equilibrium roll 
rate Pe in the models of Fig. A, 3 equal zero).wten Pe = 0, 
Qe (i.e., equilibrium pitch rate) has been formed to have 
negligible influence in the lateral model compared 

to and was therefore set equal to zero. 

B. OHCOOPLED PITCH CHANNEL AUTOPILOT 

A general block diagram of an uncoupled pitch channel 
autopilot is shown in Fig. B.1 . 

A normal acceleration command g’s) is applied to 

the pitch control law which uses measurements of missile 
body pitch angular rate (q) and pitch normal acceleration 
determine the required actuator command (Sp^) . The 
actuator is modeled as a first-order lag at 30Hz (188.4 
rad/sec) . 

The dynamic model is linearized about a trim angle-of- 
attack as described in appendix A. 

From the block diagram of Fig. B.1, the parameters of 
interest for the present analysis are given respectively 

^ C ' Cm ' • 

*^o<. 
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Control Law 1 Actuator Dynamic HooHL 







Figure B.1 Oncoupled Pitch Channel autopilot. 
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For ease notation the correspondence of state variables and 
constants are given in Table Vll 

Performing block diagram manipulations see [aef. 5], one 
yields to the following state variable equations: 

aerodynamic equations 

XI = C2 A3 X2 + C2 A4 X3 (B.1) 



X2 = XI - K Cl ai X2 - K Cl A2 X3 (B.2) 

Control law equations 

X3 = - C3 X3 + C3 Conv X6 (B.3) 



X4 



Cl C4 A1 X2 - Cl C4 X3 - C4 X4 



(B.4) 



X5 = C7 X4 - C5 X5 - C6 NZC (B.5) 

Actuator equation 

X6 = C E A7 X3 “ 6 X6 + C E as X8 + C E A6 XII (B.6) 



C. COUPLED ROLL-riW CHANNEL AUTOPILOT 

A general block diagram of a coupled roll- yaw channel 
autopilot is shown in Fig. B.2. 
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CoNFROL Law5 '-Actuators I OYNAmc Model 




Figure B.2 Coupled Roll-law Channel Autopilot. 
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The block diagram (Fig.B.2) shows, eight parameters of 
interest that are; 

A4 = , 

A5 = , A6 = , A7 = Cjj ^ , and A8 = - 

For ease of notation the correspondence of state vari- 
ables and constants is given in Table Vlll. 



Eqns.B.7 through B.17 give the state variable equations 
of the coupled roll-yaw channel autopilot: 

Aerodynamic equations 

XI = C Conv(A4 X3 + A3 X8 + A5 XII) (B.7) 



X2 = C Conv(A7 X3 + AS X8 + A6 XII) 



(B.8) 



X3 = - XI - (alpha/Conv) X2 (B.9) 

+ KB A A2 X3 + KB A A1 XII 

X4 = -8 X4 -17.6 X12 ♦ 17.6 PHC (B. 10) 

Control law equations (roll) 

9 

X5 = - (0.755/Conv) X2 - C D A7 X3 (B. 11) 

♦ (0.755 - 8 D) X4 - 5 X5 - C D A3 X8 - C D A6 XII 

- 17.6 D X12 +17.6 PHC 
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X6 



C E A7 X3 - 6 X6 ♦ C E A8 X8 «■ C E A6 XII 



(B. 12) 



X7 = - F KC (0.755/Conv) X2 - (D ♦ E) F KC A7 X 3 (B. 13) 

+ F KC (0.755 - 8 D) X4 ♦ KC(15 - 5 F) X5 

♦ KC(6F - 15) X6 - 15 X7 - (D ♦ E) F KC C A8 X3 
- (D ♦ E) F KC C A6 XII - F KC D 17.6 X12 

♦ F KC D 17.6 PHC 

* 

X8 = 188.4 Conv X7 - 188.4 X8 (B.14) 

Control law equations (yaw) 

X9 = (K1 A A2/T1) X3 - X9/T1 + (K1 A a1/Tl) XII (B.15) 



X10 = (K2/Conv) XI - (K2 H ALPHAB) /Conv X2 (B. 16) 

+ (K2/10) ((K1 A A2)/T1 + B A4 - H ALPHAB C A7 X3 

+ (K2/10) (B A3 - H ALPHAB C A8 X3 ♦ K2(1 - (1/10 T1) X9 
+ (K2/10) ((K1 A A1)/T1) ♦ B A5 - H ALPHAB C A6) XII 
Actuator equations 

XII = 188.4 conw X10 - 188.4 XII (B. 17) 



X12 = X2/COUV (B. 18) 
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D. AEEODYHAHIC DATA-LINEAR APPROXIS ATION 

A linear approach was used in the design and stability 
analysis of the aultopilots of tha pitch, yaw, and roll 
channels, both uncoupled and coupled. According, a linear 
approximation of the aerodynamic derivatives at m=3.95 was 
provided, about which the system could be perturbed. These 
linearized aerodynamic derivatives are given in [Ref- 2] and 
are presented here in Table IX. 
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TABLE ?II 

Correspondence of Syabols (Oncoapled Pitch Autopilot) 



“1 



q 



5 



p 



X 

Y 



Sp 

qs/d 

57.3qsd/I 



188.4 

150.0 

6.0 



0.53544 

0.48 

0.38375 

3.07 

0.48 

57.3 



XI 

X2 

x3 

X4 

x5 

X6 

A1 

A2 

A3 

A4 

Cl 

C2 

C3 

C4 

C5 

C6 

C7 

C8 

C9 

K 

Conv 
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T&BLE VIII 

Correspondence of Symbols (Couplsd Roll-Yaw Autopilot) 



r 

P 

fi 

X 

II 

XI 



iftc 

Y 

he 

Sy 



Cv/- 



'^SR. 



ysy 

'isa. 



XI 3.1416 S 

X2 57.3 Conv 

X3 qs/d A 

X4 qsd1/l22 B 

X5 qsdl/I^x C 

X6 0.05033 D 

X7 0.078 E 

X8 0.25 F 

X9 e ALPHA 

X10 c<e /S7.3 ALPHAS 

XII 0.48 KB 

X12 0.25 T1 

A1 0.84 K1 

A2 6.08 K2 

A3 1.0 KYP 

A4 4.17 KC 

A5 0,48 K 

A6 2.0 D1 

A7 "q QB 

A8 
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TABLE IX 

Linearized Aerodynamic Derivatives 



o< = 10° 







0. 054 







r - 


-f 


0. 024 ■ 


» 

c 

) 




Ca 




0. 027 


5(3 




c.. 





0. 015 


Y,fr 








0. 039 







c« 




0. 010 






Cw 




0. 006 








+ 


0. 014 







C - ............. 


+ 


0. 023 


Ur. 




c. 


^ 


0. 220 


Nc< 




c., 




0. 020 










0. 0137 




»ne< 




c 


__ 


0. 055 


»Ttfp 
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aPPENDII S 

NONLINEAR SYSTEM MODEL 



A. INTRODOCTION 

A general block diagram of the nonlinear model is shown 
in Fig. C.1 The three dimensional nonlinear aerodynamic model 
presented here is for the same conditions used in the linear 
model shown in appendix B. The same flight condition used 
for the linear case is the same for the nonlinear model (i. 
e. , 60 kft altitude, Mach 3,95). The control laws are the 
same used for the linear models except for a minor modifica- 
tion to the coordinating branch dependence on angle-of- 
attack and also the inclusion of anti-gravity bias as stated 
in [Ref. 2], 



B. CONTROL LARS 

The control laws used for the following nonlinear 3-D 
studies were the same as those used for the linear models, 
except for the gain o< shown in the bold line of coordina- 
tion branch in Fig. C.1. The new gain ^ is held constant 
at one degree magnitude for angle-of-attack less than one 
degree positive. For angle-of-attack greater than one 
degree positive, the gain is equal to the angle-of- 

attack. This keeps coordination for very small angle-of- 
attack. 

Gravity effects were not included in the linear models 
in appendix B, because it was assumed to have a negligible 
influence on autopilot stability and response for pertuba- 
tions about a missile trim condition. However, gravity 
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Figure C« 1 



Nonlinear Model. 
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effects were included in the following nonlinear model where 
the missile body-fixed yaw axis will be subjected to the 
full force of gravity and may, therefore, have a significant 
influence on sideslip. 

The gravity effects are compensated for by pitch and 
roll motion of th^ missile which have less influence on 
sideslip than yaw motion. In inertial rectangilar coordi- 
nates one has: 

, ^ =acceleration command in inertial direction 



= V|^^_ cos(e) 



(C.1) 



where 

*' 12 . = guidance command ( g’s ) 

-Cos© = anti-gravity bias command ( g's ) 

= acceleration command in the inertial direction 

guidance command ( g's ) 

There is no gravity effect in Y direction. In order to 
assure an anti-gravity bias of just one g, oa e needs to 
modify the anti-gravity command as follows; 

elliptical airframe inertial acceleration command 

= Nzc-0. 91 3cos (^) 

c 



C. PEBFOHANCE COHHANDS 

The commands 2 gees is first applied in the 0<> or 
upward direction at 2 seconds. since both the missile 
roll angle and roll angle command are at zero degrees, there 
is no roll motion and the missile turns upward as a 
skid-to-turn controlled missile. At 5 seconds, a second 
2 g's inertial guidance command is applied in the downward 
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or 1800 direction. The missile is commanded to roll through 
1800 while moving in a coordinated manner in yaw and roll to 
minimize sideslip angle and prevent or minimize negative 
angle -of -at tack. 

The response of achieved maneuver and cross-plane accel- 
erations during the first 2 seconds are due to initial 
conditions, gravity and anti-gravity bias effects. The 
initial conditions were added to minimize the transients 
which result when gravity bias commands the autopilot for 
constant altitude missile flight. The initial conditions 
were : 

^ = angle-of-attack = 2.41 deg 

Sjp = pitch tail angle = 0.658 deg 
© = pitch Euler angle = 3.65 deg 

output of pitch acceleration feedback lag = -1.0 g*s 
pitch control law acceleration error 
lag prior to dc gain = - 0.0105 
5pc = pitch actuator command = 0.558 deg 

The achieved maneuver plane acceleration (Eqn.C.2) is 
calculated from the body-fixed acceleration 
the roll angle {) as follows; 









^6 



and 



"I 





^ Y © S\n ^ 



(C.2) 



During the first command, achieved body-fixed yaw accel- 
eration and missile roll angle are equal to zero, 

because the roll channel is not commanded. Therefore, 
achieved maneuver plane acceleration is equal to the body- 
fixed acceleration z^. During the second command the 
missile roll angle has the same roll angle response of the 
roll autopilot. 
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D. AERODYHAMIC DATA - NONLINEAR REPRESENTATION 

The aerodynamic data are given in [Ref. 2]. The entire 
study was conducted using B=3.95 and the aerodynamic coeffi- 
cients are based on a body- fixed axis system of Pig.B. 2. 

For reference, the normal force and pitching moment 
plots given in appendix B of [Ref. 2 ] at 8=3. 95 were repro- 
duced in Table X and XI, which give just the values for 
points of interest. These tables present, too, the corre- 
spondent derivatives that were found by simulation using a 
fortran program shown in appendix F. The aerodyiamic deriv- 
ative ofCjj, , and Cy with respect to sideslip 

angle (h , yaw control 6y , and roll control are 

presented in Table XII, XIII and XIV which give just values 
for points of interest and includes the correspondent deriv- 
atives that were found by simulation using a fortran program 
given in appendix P. 

The block diagram of Fig.C.1 shows 10 parameter of 

interest that are : 

A1 = C^(oc,<fp) , A2 = , A3 = , 

A4 = C (X) , A5 = C. (X) , A6 = Cy (x) , 

A7 = » A8 = (X) , A9 = , 

and AA = Cy^ (<^) . 

Using the same procedure as given in appendix B one can 
have the Tables X7 and XVI that give the correspondence of 
symbols for the nonlinear model. Egns.C.3 through C.21 give 
the state variable equations of the nonlinear model. 
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XI = - C4 XI - C02 C4 CN(X6,X4) (C.3) 



X2 = C7 XI - C5 X2 ♦ C6 C19 C0S(X17) - C6 NZC (C.4) 



X3 = - C7 C8 X1 ♦ (C5 C8 - C9) X2 + (C9/Conv) X5 (C.5) 



+ C6 C8 NZC 

# 

X4 = - C3 X4 ♦ C3 Conv X3 (C,6) 

« 

X5 = (X17 Xl8)/Conv ♦ C1 Conv CM(X6,X4) (C.7) 

X6 = X5 - KB CM(X6,X4) - (X16 X18|/Con7 (C.8) 

♦ 

X7 = X5 Cos(Xl9) /Conv - X17 Sin (XI 9 ) /Conv (C.9) 

f 

X8 = - C10 X8 + (C12 - C11 C13) X10 (C, 10) 



XI 5 



(C. 11) 



♦ C11 C14 PHC - C11 C14 X19 - CO C1l(Cy (X6) 

♦ (x6) X13)- C12 X18/Conv ' 

X9 = - X9/T1 + (K1/T1) C02(C (X6) X16 



♦ C^^(X6) XI 5) 
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X10 



C13 X10 ♦ C14 PHC 



C14 X19 



(C. 12) 



XII = - C15 XII ♦ C16 C0(C^^(x6) X16 (C.13) 

♦ C (X6) X15 ♦ C (X6) X13) 

X12 = - C7 X12 ♦ K (C17 - C10/C18) X18 (C, 14) 



+ (K/C18) (C12 - C11 C13) X10 + K C11 (C14/C18| PHC 

- K C11 C14/C18 X19 - (K/C18) CO C (X6)(C11 

+ C16) X16 - (K/C18) CO Cjp(X6)(C11 + C16) X15 

- (K/C18) CO Cj(^(X6) (Cl 1 + C16) X13 

- (K C12/CDI1V) XI 8 + (C15/C18 - C17) XII 

# 

X13 = - C3 X13 ♦ C3 Conv X12 (C. 15) 



X14 = (K2/10) (- X9/T1 + K1 T1 C02 (C (X6) X16 (C. 16) 

+ Cy^^(X6) X15) - (X15 XI 8) /Conv2 (X5 - KB C02 CN(X6rX4) 

- (X16 X18)/Conv X18 - (KXP/Conv) (X5 - KB C02 CM(X6,X4) 

- (X16 X18)/Conv X18 - (KYP/Conv) CO XX6 (C^^(X6) X16 



- Cjj^(X6) X15 



* ^13) ♦ (K2/Conv) X17 



+ K2 X9 - (K2 KYP/Conv2) XX6 X18 



X15 = - C3 X15 ♦ C3 Conv X14 



(C. 17) 



X16 = KB 



C02 (Cy^ (X6) 



X16 * 



(X6) 



X15 



+ (X16 X18) /Conv - X17 



(C. 18) 
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X17 = - (X5 XI 8) /Conv ♦ C01 Conv(3^^(X6) X16 
+ (X6) X13) 

X18 = CO Conv (Cn (X6) X16 + C„ (X6» X15 

* C (X6) X13) 



XI 9 = X18/Conv 



(C. 19) 



(C. 20) 



(C. 21) 
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T&BLE X 

Normal Force Coefficient and Derivatives M=3.95 









^ Cfj 


^ C )S4 <fp) 

1) <£|^ 


-4.00C0 


-lO.OOOC 


-1.C500 


0.1842 


0.0308 


-4. COCO 


0.0000 


-0.8000 


0 .1453 


0 .0208 


-4.C0CC 


5.000C 


-0.7000 


0.2496 


0.0196 


-4.0000 


10.0000 


-0.6000 


0.2141 


0.0208 


0.0000 


-10.0000 


-0.2500 


0 .2110 


-0.0332 


0.0300 


o.oooc 


-o.caoo 


0.20 74 


0.0448 


0.0000 


5.0000 


0.1300 


0.1827 


0.0336 


0.0000 


lO.OOOC 


0.2000 


0.1929 


-0.0112 


4.0000 


-10.0000 


0.6000 


0.2092 


0.0200 


4.0000 


0.0000 


0.8000 


0 .2253 


0.0200 


4. COCO 


5.000C 


c.sooo 


0.2196 


0.0200 


4.0000 


lO.OOOC 


l.COOO 


0.2141 


0.0200 


8. COCO 


-lO.OOOC 


1 .4500 


0 .2272 


0 .0042 


8.0000 


O.OOOC 


1.7000 


0.2265 


0.0342 


8.0000 


5.0000 


1.6500 


0.2288 


0.0229 


8.0000 


lO.OOOC 


l.SOOO 


0.2256 


-0.0058 


• 12.0000 


-lO.OOOC 


2.4000 


0.2321 


0.0200 


12.0000 


0.0000 


2.6000 


0 .2188 


0.0200 


12. COCO 


5.000C 


2.7000 


0.2153 


0.0200 


12.0000 


lO.OOOC 


2.6000 


0.2335 


0.0200 


16.C0CC 


-lO.OOOC 


3 .3000 


0 .2319 


0 .0083 


16.0000 


0. OOOC 


3.50C0 


0.2482 


0.0283 


16.0000 


5.0000 


3.6500 


0.2599 


0.0308 


16.0000 


lO.OOOC 


3.6000 


0.2654 


0.0283 


20.0000 


-10.0000 


4.3000 


0.2653 


0.0300 


20. COCO 


o.oooc 


4.6000 


0.2884 


0.0300 


20.0000 


5.0000 


4.7500 


0.2825 


0.0300 


20.0000 


10.0000 


4.5000 


0 .2798 


0.0300 


24.0000 


-lO.OOOC 


5.40CO 


0.2819 


0.0300 


24.0000 


0.0000 


5.7000 


0.2482 


0.0300 


24.00CC 


5.0000 


5.6500 


0 .2599 


0.0300 


24.0000 


lO.OOOC 


6.C000 


0.2654 


0.0300 



I 



J 
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TiBLB XI 

Pitching Moment Coefficient and Derivatives M=3.95 






o< 


S 




ip) 

lioL 




-4.000C 


-lO.OOOC 


0.4850 


0.0130 


-0.0159 


-4.0000 


-5.0000 


0.2800 


-0.0200 


-0.C61C 


-4.0000 


0.0000 


-0.C750 


0 .0190 


-0.0759 


-4.00C0 


lO.OOOC 


-0.6300 


0.C197 


-0.0149 


0.0000 


-10.0000 


0.5500 


0.0185 


-0.055C 


o.oocc 


-5.000C 


0.2750 


0 .0131 


-0.0550 


0.0000 


0. 0000 


-o.cooo 


0.0186 


-0. 0550 


0.0000 


10.0000 


-0.5500 


0 .0195 


-0.0550 


4.00CC 


-lO.OOOC 


0.6250 


0.0180 


-0.0551 


4.0000 


-5.000C 


0.2500 


0.0200 


-0.055C 


4.0000 


0.0000 


0.C750 


0 .0190 


-0.0551 


4.0000 


lO.OOOC 


-0.4800 


0.0147 


-0.0561 


8.000C 


-10.0000 


0.7100 


0.0294 


-0.0563 


8. 00 CO 


-5.000C 


0.^300 


0.0232 


-0.0558 


8.0000 


O.OOOC 


0. 1500 


0.0178 


-0.0563 


8.0000 


10.0000 


-0.4300 


0.0118 


-0.0603 


12.000C 


-lO.OOOC 


0.6500 


0.0331 


-0.0640 


12.0000 


-5.0000 


0.5250 


0.0184 


-0.0655 


12.0000 


0.0000 


0.2000 


0 .0034 


-0.0640 


12.0000 


10.0000 


-0.2900 


0.0057 


-0.052C 


16.C000 


-10.0000 


0.9500 


0 .0182 


-0.0725 


16.0000 


-5.000C 


0.5750 


0.0120 


-0.0763 


16.0000 


0.0000 


0.2000 


0.0059 


-0.0725 


16.00CC 


lO.OOOC 


-0.4000 


-0 .0122 


-0.0425 


20.0000 


-lO.OOOC 


l.CCOO 


0.0065 


-0.0756 


20.0000 


-5.0000 


0.6220 


0 .0065 


-0.0755 


20.00CC 


O.OOOC 


0.2450 


0.0067 


-0.0752 


20.0000 


10.0000 


-0. 5000 


-0.0395 


-0.0736 


24.00C0 


-10.0000 


1 .COOO 


-0 .006 8 


-0.0841 


24.0000 


-5.000C 


0.5870 


-0.0290 


-0.0818 


24.0000 


0.0000 


0.1750 


-0.0516 


-0.0837 


24.00CC 


lO.OOOC 


-0.7300 


-0.0772 


-0.1001 
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TABLE XII 

Sideslip Derivatives and Derivatives 



«a< 


Cnt'-) 

V 


'"t) oc 


-4.0000 


0.0110 


-0.0027 


0.0000 


0.0000 


-0.00 27 


4. 0000 


- 0 . Clio 


-C. 0022 


5.2000 


-0.0135 


-0.0024 


8.0000 


-0. 0220 


-0.0027 


10.0000 


-0. 0270 


-0.0021 


12.6000 


-0.0310 


-0.0015 


16.0000 


-0. 0360 


-0.0 014 


20.0000 


-0.0410 


-0.0011 












C< 



-4.00C0 
0.0000 
4.0000 
5.20CC 
8.00C0 
10.0000 
12.60CC 
16.0000 
20. COCO 



0.0240 
0. 0240 
0. 0240 
0.0240 
0. 0240 
0.C240 
0. C270 
0. 0290 
0.0320 



C.OOOO 
C.0000 
G.OOCO 
0.0000 
C.00C0 
C.0005 
0.0 009 
C. 0 007 
0.0008 



0< 



Cy {.e^ ) 



^ Cy 






c>4 



-4.0000 

0.00 00 

4.0000 
5. 2000 

8.0000 

10.0000 
12.6000 
16.0000 

20.0000 



•0. 0490 
•0. 0440 
-0.0480 
•0. C490 
•0.0520 
•0. 0550 
■0. 0580 
•0. 0610 
•0. 0650 



0.0 024 
0.0001 
-0.0009 
-0.0 009 
-0.0013 
-0.0013 
- 0.0010 
-0.0009 
- 0.0011 
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TABLE XIII 

Yaw Control Derivatives aod Derivatives 










D wL. 


-4.0000 


0. 0025 


-C. 00C7 


0.0000 


0.0000 


-0.0006 


4.00CC 


-0. C020 


-0.0 014 


5. 2000 


-0.0040 


-0.0015 


8.00C0 


-0.0075 


-0.0013 


10. COCO 


-0. ClOO 


-C.0 014 


12.6000 


-0.0140 


-0.0015 


16.C000 


-0. C190 


-0.0014 


20.0000 


-0. C240 


-0.0011 


o< 


Ch 


^ oc 


-4.00C0 


-0. 0400 


0.0000 


0.0000 


-0. 0400 


C.0 0C0 


4.0000 


-0.0400 


C.OOOO 


5.20CC 


-0. C400 


0.0000 


8.0000 


-0.0400 


C.00C0 


lO.OOCO 


-0. 0400 


C.OOOO 


12.600C 


-0. C400 


-C.00C3 


16.0000 


-0.0420 


-0.0007 


20.000C 


-0. C453 


-0.0010 



oC 















-4.00CC 

0.0000 

4. 0000 
5.20C0 

8.0000 
10. COCO 
12.60CC 
16.0000 
20.0000 



0. C150 
0.0150 
0.0150 
0. C150 
0.C150 
0. C150 
0. 0150 
0.C170 
0. C195 



0.0000 
C.00C0 
0.0000 
C.OOOO 
C.00C0 
0.0000 
C.0 003 
C.0006 
0.0 006 
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TABLE XIV 






Roll Control 


Derivatives and 


Derivatives 
















c) 




- 4.0000 


0.0232 


0.0000 




O.COCO 


0. 0232 


0.0 0 00 




4.0000 


0. 0232 


0.0000 




5. 2000 


0. 0232 


C . 0 000 




8.0000 


0.0232 


0.0 0 00 




10.0000 


0 . 0232 


0.0 003 




12.6000 


0. 0250 


0.0006 




16.0000 


0. 0270 


0.0006 




20.0000 


0 . 0298 


C . 0 008 












o< 


op. 


oc 




- 4 . 00 CC 


- 0.0045 


0.0011 




o.oocc 


0 . 0000 


0 . 0 011 




4.0000 


0.0045 


0.0 009 




5. 2000 


0. 0055 


0.0012 




8.0000 


0. Clio 


C . 0 018 




1 0.00 00 


0.0145 


0.0 018 




12.6000 


0. 0195 


0 . 0 018 




16.0000 


0.0250 


0.0017 




20.0000 


0. 0320 


0.0018 
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TABLE XV 

Correspondence of Symbols ( nonlinear Model ) 
XI p 




T 


X2 


‘f 






X3 









X4 







q 


X5 







of. 


X6 






b 


X7 


1 Sy 






T1 


X8 







Y2 


X9 


'K 




XI 


X10 


> 

^ 




X2 


XII 




Si^c 


X12 







Sr. 


X13 


qSd/Ip^^ •••••••«• 






X14 


qSd/Iz2 






X15 


qSd/W 




P> 


X16 






r 


X17 


3.1416 


J 
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TABLE XVI 

Correspondence of Symbols (Coat. Nonlinear Sodel) 



188.4 C3 

150.0 C4 

6.0 C5 

0.530544 C6 

0.48 C7 

0.38375 C8 

3.07 C9 

5.0 CIO 

0.05033 C11 

0.755 C12 

8.0 C13 

17.6 C14 

6.0 C15 

0.078 C16 

15.0 C17 

4.0 C18 

0.913 C19 

0.25 T1 

0.839 K1 

6.08 K2 

1.0 KYP 
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UNCOUPLED PITCH AUTOPILOT - CSMP PROGRAM 
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AERODYNAMIC EQUATIONS 
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CONTROL LAW EQUATIONS 
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